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Abstract 

We calculate the electro- and chromomagnetic dipole coefficients and in 

the context of the minimal Randall-Snndrum (RS) model with a Higgs sector localized 
on the IR brane nsing the five-dimensional (5D) approach, where the coefficients are 
expressed in terms of integrals over 5D propagators. Since we keep the full dependence 
on the Yukawa matrices, the integral expressions are formally valid to all orders in 
In addition we relate our results to the expressions obtained in the Kaluza- 
Klein (KK) decomposed theory and show the consistency in both pictures analytically 
and numerically, which presents a non-trivial cross-check. In Feynman-’t Hooft gauge, 
the dominant corrections from virtual KK modes arise from the scalar parts of the W^- 
boson penguin diagrams, including the contributions from the scalar component of the 
5D gauge-boson field and from the charged Goldstone bosons in the Higgs sector. The 
size of the KK corrections depends on the parameter y*, which sets the upper bound 
for the anarchic 5D Yukawa matrices. We find that for y* > 1 the KK corrections are 
proportional to yj. We discuss the phenomenological implications of our results for the 
branching ratio Br(i? —>■ the time-dependent CP asymmetry Sk*'^, the direct CP 

asymmetry A^^'^ and the CP asymmetry difference We can derive a lower 

bound on the first KK gluon resonance of 3.8 TeV for y* = 3, requiring that at least 
10% of the RS parameter space covers the experimental 2cj error margins. We further 
discuss the branching ratio Br(.B — XsDl~) and compare our predictions for C'y.^^gqo 
and with phenomenological results derived from model-independent analyses. 


1 Introduction 


In July 2012 the Higgs boson, the last missing piece of the Standard Model (SM), was discov¬ 
ered at the Large Hadron Collider (LHC) at CERN [1, 2]. Since then the hierarchy problem, 
i.e. the question about the mechanism that stabilizes the Higgs mass near the electroweak 
scale, is no longer a hypothetical issue. A promising possibility to solve the hierarchy problem 
is offered by Randall-Sundrum (RS) models [3], in which the SM is embedded in a slice of 
anti-de Sitter space while the Higgs sector is localized on the “infra-red (IR) brane”, one of 
two sub-manifolds bounding the extra dimension. The smallness of the electroweak scale can 
then be explained by the fundamental ultra-violet (UV) cutoff given by the warped Planck 
scale, whose value near the IR brane lies in the TeV range. Moreover, by allowing the fermion 
helds to propagate in the bulk, these models provide a natural explanation for the hierar¬ 
chies observed in the flavor sector [4-6] and the smallness of flavor-changing neutral currents 
(FCNCs) [7-13]. 

In this paper we investigate the FCNC process 6 ^ sy in the minimal RS model with 
a brane-localized Higgs sector. For two reasons this transition is very interesting in order 
to search for new physics. In the SM the dipole coefficients are one-loop suppressed and 
the transition is logarithmically suppressed by the GIM mechanism [14]. In order to include 
the effects of the RS model on the transition 6 —)■ sy we implement an effective Lagrangian, 
in which the heavy Kaluza-Klein (KK) quarks and bosons are integrated out. The most 
important operators are the electromagnetic dipole operators 

Qi, = Qj, = -^rsa.^F^^PLb, ( 1 ) 


with = I ['Jfj.j'yu] and the projection operators Pr^l = |(1 ± ys)- Due to operator mixing 
we also consider the chromomagnetic dipole operators 
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S G'^^'taPL b , 


( 2 ) 


where ta are the generators of SU{3)c- The main focus of our paper lies on the derivation 
of integral expressions for the dipole coefficients at the one-loop level using hve-dimensional 
(5D) propagators in the mixed position-momentum space and with the full dependence on the 
Yukawa interactions imposed by the mixed boundary condition at the IR brane. 

In the literature, the first discussions on 6 —)■ sy in the RS model can be found in [7, 8, 15]. 
There, the authors claimed that the penguin diagrams with the exchange of charged Higgs 
scalars (Goldstone bosons of the boson) along with KK fermions gives the dominant 
contribution to the dipole coefficients. The diagram with the exchange of KK gluons was found 
to be approximately aligned with the 4D down-type Yukawa matrix and therefore subleading. 
Furthermore, the authors claimed that the dipole coefficients in the brane-localized Higgs 
scenario were logarithmically divergent and sensitive to the UV cutoff. It was shown in [16] 
that the diagrams contributing to the leptonic decay /x —)■ ey at one-loop are indeed hnite. 
With the same technique the authors of [17] investigated the process 6 —)■ sy working with 5D 
propagators and treating the Yukawa interactions as perturbations. In [18], one of us discussed 
the 6 —)■ sy process in the minimal RS model with a brane-localized Higgs sector working in 
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the KK-decomposed theory, where the dipole coefficients are expressed via inhnite sums over 
the contributions from different levels of KK modes. In [19] the authors calculated the dipole 
coefficients in the custodial RS model with a brane-localized Higgs sector, focusing only on 
the diagrams with an exchange of the first level of KK fermions along with gluons and charged 
Goldstone bosons. Recently, the authors of [20] studied lepton flavor violation in RS models 
in the 5D framework, where they discussed the electromagnetic (leptonic) dipole operator in 
RS models with a brane-localized or nearly brane-localized Higgs and treated the Yukawa 
interactions as perturbations. In the present work, we perform a complete calculation of the 
electro- and chromomagnetic (quark) dipole coefficients including all contributions at one-loop 
order in the minimal RS model with a brane-localized Higgs sector. We derive expressions for 
the dipole coefficients using 5D propagators computed by retaining the full dependence on 
the Yukawa interactions. In contrast to [17, 20], we derive 5D expressions that are formally 
valid to all orders in In contrast to [20], we focus on the quark dipole coefficients, 

including the contributions of the chromomagnetic dipole operator. In addition, we derive 
formulas in the KK-decomposed (4D) theory including the contributions from all KK levels 
and show the consistency with the results obtained in the 5D framework. 

After introducing the model and setting up the notation in Section 2 we derive formulas for 
the dipole coefficients in the 5D framework in Section 3. In Section 4 we compare our results 
with the expressions in the KK-decomposed theory and analyze the different KK contributions 
to the dipole coefficients. After implementing the renormalization-group (RG) evolution from 
the KK scale down to the R-meson scale we discuss the phenomenological implications in 
Section 5. Our main results are summarized in the conclusions. 


2 Theoretical setup 

We focus on RS models where the electroweak symmetry-breaking sector is localized on or 
near the IR brane. The extra dimension is chosen to be an /Z 2 orbifold parametrized by 
a coordinate 0 G [—7r,7r], with two 3-branes localized on the orbifold fixed-points 0 = 0 (UV 
brane) and |0| = tt (IR brane). The RS metric reads [3] 

ds^ = 7]^^ dx^dx'^ — r'^dcj)^ = ^ (ri^„ dx^dx'^ — dt^^ , (3) 

t \ ^KK ) 

where with ct( 0) = /cr|0|, is referred to as the warp factor. The size r and curvature 

k of the extra dimension are assumed to be of Planck size, k rsj \ IT" Mpi. The quantity 
L = (T{7r) = km measures the size of the extra dimension and is chosen to be L ~ 33 — 34 
in order to explain the hierarchy between the Planck scale Mpi and the TeV scale. We dehne 
the KK scale Mkk = ke^ with e = which sets the mass scale for the low-lying KK 

excitations of the SM particles. On the right-hand side of (3) we have introduced a new 
coordinate t = e whose values on the UV and IR branes are e and 1, respectively.^ In our 
analysis we consider the minimal RS model, adopting the conventions and notations of [10]. 

^The dimensionless variable t is related to the conformal coordinate z frequently used in the literature by 
the simple rescaling z = I/Mkk = R't- 
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The gauge group is taken to be SU{3)c x SU{2 )l x U{1)y like in the SM, and it is broken to 
SU{3)c X U(l)em by the Higgs vacuum expectation value (vev). 

A tree-level analysis of electroweak precision observables, mainly the T parameter, implies 
that the mass of the hrst KK gluon resonance is pushed to values Mg(i) > 11.3 TeV at 95% 
conhdence level (CL) [21], where we have used the most up-to-date values from [22], Loop 
corrections could potentially change this bound in a signihcant way. Also, it is conceivable 
that new-physics contributions arising in a UV completion of the RS model bring S and 
T back into the phenomenologically favored region. Nevertheless, one usually considers RS 
models with a built-in protection for the T parameter by implementing a custodial symmetry 
via the gauge group SU (3)c x SU{2)i x SU (2)/j x U (l)x x [23-25]. Then the bound from 
electroweak precision observables reduces to M^(i) > 4.8 TeV at 95% CL [29]. But on the other 
hand, the sensitivity of Higgs physics on virtual effects from heavy KK excitations is strongly 
increased in the model with custodial symmetry, due to the enlarged fermion multiplicity for 
each KK level. Comparing predictions for the signal rates of the Higgs decaying into pairs of 
electroweak gauge bosons with data from the LHC excludes KK gluon resonances lighter than 
(15 — 20) TeV X (|/*/3) at 95% CL [27], where the precise value depends on the details of the 
localization of the Higgs sector near the IR brane. Here sets the upper bound for the entries 
of the anarchic 5D Yukawa matrices, |(Y)j)jj| < y^. In the minimal RS model the resulting 
bounds are much weaker. For values of y^ > 0(1), the custodial RS model thus loses its main 
advantage of allowing for lighter KK resonances such that the minimal RS model is just as 
promising nowadays. Note also that the parameter measuring CP violation in kaon mixing 
requires KK gluon masses in the range of 10 TeV (with moderate hne tuning) irrespective of 
whether the minimal or the custodial RS models are considered [9, 26]. 

Higgs localization 

In this work we focus on the RS model with a brane-localized Higgs field, where the inverse 
characteristic width A/^ of the Higgs held along the extra dimension is assumed to be much 
larger than the inherent UV cutoff near the IR brane, i.e. A/j ^ Atbv ~ several Mkk [28]. 
It is well known that quantum helds can be strictly localized on orbifold hxed points, and 
in such a scenario the quantity A^ can indeed be inhnite or arbitrarily large. Similar to the 
case of Higgs production via gluon fusion [29, 35-45], we will hnd that our results for the 
Higgs contributions to the 6 —)■ sy and b ^ sg dipole coefficients are sensitive to details of 
the localization mechanism. For these contributions we sometimes extend our results to the 
case of a so-called narrow bulk-Higgs scenario, where the Higgs held lives in the bulk with an 
inverse width such that Mkk ^ A^ Axev [29]. This is a special case of a general bulk-Higgs 
with Ah ~ V, where v is the Higgs vev. 

The authors of [30, 31] have calculated the RS contributions to the dipole operators in the 
KK-decomposed theory for a general bulk-Higgs held, where the localization parameter is taken 
to be /9 ~ 1 (in our notation (3 ~ Ah/v). Numerically, they also discuss the quasi IR-localized 
limit by increasing [3, i.e. by pushing the Higgs prohle towards the IR brane. They hnd that 
heavy KK fermion modes with masses fUq^ ~ /3Mkk yield unsuppressed contributions in the 
case where the Higgs inverse width is of order the UV cutoh, (3 ATev/'i' {^h ~ Axev)- In 
this case, there are still some high-momentum KK modes that can probe the “bulky nature” 
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of the Higgs field. In addition, the authors of [20, 30] have observed the non-decoupling of 
heavy KK excitations of the Higgs boson itself in the quasi IR-localized limit of large f3. These 
hndings show that the results of Higgs-induced contributions to the dipole operators depend 
on the implementation of the Higgs sector (see [20] for a nice recent discussion of this point). 
In our brane-localized Higgs scenario, where A/j 3> AxeV) we observe that the 5D description 
of the dipole coefficients gives results which are consistent with the description in the KK- 
decomposed theory when summing over the hrst few KK levels. This shows that heavy KK 
modes near the UV cutoff decouple. Furthermore, KK excitations of the Higgs doublet do not 
arise in this version of the model. 


Gauge sector 

In the minimal RS model the SM gauge group lives in the bulk and is broken to U (1)em on the 
IR brane, where the Higgs held develops a vev. Details for the implementation of the Higgs, 
gauge-boson, and gauge-hxing sectors in the context of this model (and using our notations) 
have been given in [10]. The KK-decomposition for the 5D gauge-boson held BM{x,t) is in 
general given by (with B = A,G, W, Z) 

B^ix^t) = dtXnit), (4) 

where Bj^^ are the KK modes of the gauge bosons with masses m^. The scalar particles 
are “unphysical” in the sense that they provide the longitudinal degrees of freedom 
of the W, Z bosons (n = 0) and their KK modes (n > 1), and thus they can be gauged away. 
Similarly the scalar particles provide the longitudinal degrees of freedom for the 

photon and gluon KK modes. The scalar helds and Z^ mix with the charged Goldstone 
bosons arising from the Higgs sector. Assuming for the time being that the scalar sector is 
localized on the IR brane, we parameterize the Higgs doublet after electroweak symmetry 
breaking in the usual form 


y/2 \v + h{x) + i(p^{x) 


where v denotes the Higgs vev in the RS model. We determine the vev v from the shift to the 
Fermi constant Gp, which can be derived in the RS model by considering (at tree level) the 
effect of the exchange of the inhnite tower of KK W bosons on the rate for muon decay. We 
hud that V = (\/2Gf)~^^'^ 1 + + O [29] with the experimentally measured value 

(\/2Ge)~^'^^ ~ 246.2 GeV. The decomposition of the scalar helds 99 ^, in (5) into the mass 
eigenstates reads [10] 
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where Cy^ = cosOyj is the cosine of the weak mixing angle, and rhwi'i^z are the leading 
contributions to the W^- and Z-boson masses in an expansion in powers of such 

that mw,z = mw,z 1 - ^ ^ (a^) SU{2)l and hypercharge 

5D gauge couplings are denoted by 5^5 and g'^. In the context of RS models the weak mixing 
angle can be expressed as = sin^ 6 *^ = g'^l{g\ + g'^')^ which can be studied experimentally 
via the Z-pole polarization asymmetries observed at LEP. We take and mw as input 
values which implies that mz is a derived quantity given in the RS model by mz{mw, -s^) = 

1 — ^ (L — 1 + + O ■ Since the prohle of the zero mode is flat up to 

corrections of order it follows that in ( 6 ) is close to 1 , and hence the 

fields 99 ^, (p^ coincide with 99 ® to leading order. We mention that one can adjust the 

gauge-hxing Lagrangian so as to cancel any mixings between the vector and scalar helds [10]. 

Quark sector 

In the quark sector the minimal RS model contains an SU{2)l doublet held Q{x,t) and 
two SU{2)l singlet helds u{x,t) and d{x,t) in the 5D Lagrangian, each of which are three- 
component vectors in generation space. The 5D fermion states can be described by four- 
component Dirac spinors [4, 5]. We use a compact notation, where we collect the left- and 
right-handed components of the up- and down-type states into six-component vectors Ua = 
{Ua,ua)'^ and Va = {DA,dA)'^ with A = L,R, which are collectively referred to as Ql,r- 
Their decomposition into 4D KK modes reads 

QA{x,t) = QA\'t)(!A\^)'^ A = L,R. (7) 

n 

The superscript n labels the diherent mass eigenstates in the 4D ehective theory, such that 
n = 1,2,3 refer to the SM quarks, while n = 4,... ,9 label the six fermion modes of the 
hrst KK level, and so on. The functions denote the wave functions of the left- and 

right-handed components of the KK mass eigenstate along the extra dimension. The upper 
(lower) components of include the prohles of the SU{2)l doublet (singlet) quark helds. 

3 Calculation of the dipole coefficients 

Like in the SM, the leading-order contributions to the 6 —)■ sy and b ^ sg dipole coefficients 
in the RS model are loop suppressed, because there are no havor-changing couplings that can 
induce a chirality hip. But in contrast to the SM there are more one-loop diagrams to be 
considered. Besides the additional exchange of KK bosons, new topologies appear due to 
the havor-changing couplings of the Higgs boson, the Z boson and its KK modes, and the 
photon and gluon KK modes. Figure 1 shows all relevant Feynman diagrams contributing 
in a general R^ gauge. Internal scalar lines of the diagrams (II), (III) and (IV) include the 
contributions from the scalar component of the 5D gauge bosons and from the corresponding 
Goldstone bosons in the Higgs sector. In this section the Wilson coefficients and C 7 'y, 8 g 


5 





Figure 1; Diagrams contributing in the minimal RS model to the transitions b ^ s'y and 
b ^ sg at the one-loop level. Solid lines denote the exchange of up- or down-type quarks 
while wavy or curled lines denote the exchange of (vector) gauge-bosons. Apart from 
diagram (I) a scalar (dashed) line includes the contribution from the hfth component of 
the gauge boson and the corresponding contribution from the Goldstone bosons in the 
Higgs sector. The extra-dimensional coordinates of the vertices are labelled according 
to diagram (I). 


are dehned via the general parametrization of the transition amplitude 




C'7'y^Sg {s'y\Q'i^^gg\b') -\- Cj^^Sg {^'y\Q7^,8g\b) 


( 8 ) 


where Gp is the Fermi constant, and Xt = V^lVtb is the relevant product of entries of the CKM 
matrix. The matrix elements in (8) are given by (Qr^) = (emfe/27r^) e*^{q) u{ps) icr^u Q^Pr u{pb) 
and {Qsg) = {gs^b/^Ti'^) e*^{q) u{ps) iCg^uQ^PRU^Pb), where the outgoing photon (gluon) mo¬ 
mentum is q = pb — Ps- The chirality-flipped matrix elements {Qr^^sg) are given by analogous 
expressions with Pp —)■ Pp. Working in Feynman-’t Hooft gauge (.^ = 1), we compute each 
amplitude in Figure 1 using the Feynman rules of the 5D theory collected in Appendix A. 

As an example we consider the penguin diagram (Ila) in Figure 1, in which a 5D W^- 
boson propagator and two 5D quark propagators arise. The corresponding amplitude with an 
external photon is given by 


y.PwrS‘‘(tp'vp.-kpi‘‘S'‘(t\t'-,p,-k)-,dpJvf(f)PL+vf(t')P,\u(tJi), (9) 


where Qu = 2/3 is the electric charge of the exchanged up-type quarks in the loop. The 
functions and V^^(t) with A = L,R denote the prohles of the physical strange- and 
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bottom-quark mass eigenstates, respectively, as defined in (7). In the above equation 65 is 
the 5D electromagnetic coupling, while represents the 5D SU{2)l gauge coupling. The 4D 
electromagnetic coupling can be obtained by e = e5/\/27rr. The 2x2 matrix Pw = P+ = 
diag(l,0) originates from the 5D Feynman rule for the vertices (with A = L,R) in 

(A.l) and projects out the profiles of the SU(2)l doublet quark fields. The 5D IF=*^-boson 
propagator in (9) can be decomposed as 


t'] k) = Bw{t, t'- -k‘^ - fO) + Bw{t, t'] -k^/i - fO) 


k^ky 

k^ 


( 10 ) 


with the KK representation 


Bwit, t'; —k"^) — ^ 

n 


- /i;2 - fO ■ 


( 11 ) 


The propagator function (11) can be calculated in closed form, see [28] for more details on the 
derivation and the solution. The 5D quark propagators in (9) can be decomposed into four 
functions differing in chirality and Lorentz structure [16, 21, 32-34], 


[AU(«, ^ +A®,Ji.i';-«;=)] + (L « fl), 


( 12 ) 


where q = u,d. The KK representations of the propagator functions read 

AiaM'; -k^) = E Sl’(‘) 2k*(0. 

n In 

-k^) = E a'S\t) Si”>'(0 . 

n 


(13) 


and analogously for and Each propagator function is a 6 x 6 matrix. The subscripts 

denote the handedness of the incoming and outgoing quark fields, such that the propagator 
function implies a chirality flip. Explicit expressions for these functions are given in 

Appendix B for the brane-localized Higgs scenario and the case of a narrow bulk-Higgs. 

Next we outline some of the basic steps needed to extract the dipole coefficients from the 
diagrams in Figure 1: 


• We perform a Taylor expansion of each 5D propagator in the external momenta Ps,Pb 
and keep the terms up to second order, since higher orders would contribute to higher¬ 
dimensional operators and yield suppressed contributions. For instance, for a 5D quark 
propagator function we apply the expansion (g = u, d) 


A%{t,t';-{p,-kr) 


‘^{Pi ■ ^ 2 





-e), 

(14) 


where Pi = Ps^b, k is the loop momentum, and A,Be {L,R}. We need to expand up 
to second order in the external momenta in order to obtain the leading effects of the 
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dipole Wilson coefficients, since the matrix elements of the dipole operators contain 
the bottom mass mb and the momentum difference q = pt — Ps- In fact, the term 
linear in pi in (14) contributes only in the RS model, and not in the SM, to the dipole 
Wilson coefficients. Analogously we can expand the 5D vector-boson propagator function 
Bsit, t'; —{pi — kY) with subscript B = A, G, hh, Z. 


• The extra-dimensional integration of the vertex with the external photon or gluon can 
be performed analytically by using the flatness of their prohles Xo '^ = 1 /and the 
orthonormality conditions for fermion and boson prohles 


dt Qfit) = 


(^Ot/ 


Stt f diy R / \ B / \ f 
nn^ -^6x6 5 Xn \^) Xn'\^) 


(15) 


where Q = Vl,V and A = L,R in the left equation and B = A,G,W, Z on the right. 


• The previous two bullets allow us to combine two 5D propagators of the same type 
when we expand them in the external momenta Ps,Pb and perform the extra-dimensional 
integration of the vertex that couples to the external photon or gluon. For instance we 
can apply (g = u, d) 


df AUt, t";-iPs - kf) F; -{Pb - kf) 


d , , ,fd 

+ k-[Ps + Pb) 


dk"^ 
2 


dk"^ 


(16) 


k ■ Psk ■ pb + {ps- ky + {pb ■ ky 


d 

dk'^ 




where we neglect terms of order (k ■ Ps)'^{k ■ Pb)^' with n -|- n' > 3. Analogous relations 
can be derived for products of different fermion and boson propagator functions. Equa¬ 
tion (16) can be used to reduce each amplitude by one extra-dimensional integration 
and one 5D propagator. 

• We perform a Wick rotation to Euclidean momenta with k^ = ik^ and kE = y—W. 

• For the matching procedure on the dipole operators we hrst use that the photon or gluon 
is on-shell, q^e^*{q) = 0, which allows us to rewrite pky*Yq) = \{p>b + Ps)^e*(g). Then, 
we can use the Dirac equation fyi{pb) = mbu{pb) and apply the Gordon identity 


U{ps) quPL,Ru{pb) = U{ps) {Ps + Pb)^PL,R - 7^ {msPL,R + mbPR,L) u{pb ), (17) 


in order to extract the Wilson coefficients. 


In the following three subsections we discuss the gauge-invariant subsets of the diagrams shown 
in Figure 1. 



3.1 Higgs contribution 

We begin with the first diagram (I) in Figure 1, in which the Higgs boson and two down-type 
5D quark propagators are exchanged. The Yukawa interactions of the Higgs boson with two 
down-type quarks are given by 






h{x) 

71 “ 


[Dl{x, t) Yd dnix, t) -h Dr{x, t) Yd dL{x, t) + h.c.] , 


(18) 


where the hrst term is often referred to as the “correct-chirality Higgs coupling” in the lit¬ 
erature, since it is also present in the SM. On the other hand, the second term couples a 
right-handed SU{2)l doublet quark held to a left-handed SU{2)l singlet, which is not allowed 
in the SM and is thus called the “wrong-chirality Higgs coupling”. The function — 1) 
denotes the normalized Higgs prohle along the extra dimension, which we take to be the reg¬ 
ularized (5-function. For the calculations we use a square box of width r] and height I/?], such 
that 

+ with (19) 

V Mkk 

where t] is related to the inverse Higgs width by t]/v. The brane-localized Higgs scenario 

corresponds to values of rj y^n/ATev, while the narrow bulk-Higgs scenario implies values 
in the range y^v/ATeY rj y^v/MKK- Note that the shape of the regularized profile is 
irrelevant as long as ?7 1. 

With the Feynman rules in Appendix A and the basic steps outlined in the beginning of 
this section we can derive an expression for the Wilson coefficient and hnd 


'^ 77,83 


K 


77-83 


dki 


V Jq kl + ml 


u3 

'^E a ^E a2 

O o 


rpd 

^RL 


{kj 


mb 


U2 u3 

'^E o _ ^ 02 '^E o3 

32 32 96 


'^Rni^E) 


Ml 


KK 


( 20 ) 


where mt is the Higgs mass. Concerning the derivatives we use the notation dkE = d/dkE- Due 
to the parametrization of the amplitude in (8) we have to divide the dipole coefficient by At 
and Gp. The couplings are given by = Qd and = 1, where Qd = —1/3 is the electric 
charge of the exchanged down-type quarks. The dimensionless propagator functions in (20) 
are dehned via 


r))i(li) = ^ dtdt' - 1 ) V™{t) Ml i'; kl) Ml cg>(i'), 


^7(^1) — 


-vMkk /■' 

72 I 


( 21 ) 


dtdt'S’’(t-l)S’’(t'-l)V™(,t)Ml AU(t,t';kl)Ml*vf(.t'), 


including the regularized (5-functions (19) and the matrix = P 12 Yd + P 21 Y^d ■ 

jector Pij for i, j = 1, 2 is a 2 x 2 matrix with zero entries except for the ij-component, which 
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equals 1. In order to perform the integrations over t and t' we need the solutions for the 
external quark prohles and the 5D quark propagators in the region t, t' E [1 — 7], 1], 

The presence of the (5-function regulator (19) implies that in the region near the IR brane, 
for t E [1 — ? 7 , 1 ], the quark prohles are determined by the coupled differential equations (with 
Q = U,T> and q = u,d) [35] 

U + jM‘ + ^ mA Q‘S\t) = , 

^ ^ ( 22 ) 

(a. - i ^xj-) sS”’(() = , 


with Xq^ = niq^/M-^Yi and g = v/{\/2 Mxk)- These equations can be simplihed for the consid¬ 
ered limit f] v/M-^k- First, the term = P+ cq — P_ Cq, which contains the bulk mass 
parameters cq and Cg, is parametrically suppressed for rjCQ^^di ^ y*d and can therefore be ne¬ 
glected. The projector P_ projects on the lower components and is given by P_ = diag(0,1). 
Secondly, the mass-dependent terms on the right side of (22) are suppressed for the SM 
quarks, since rjXg^ y^g. With these approximations, the basic solutions are given by the 
trigonometric functions (for q = u,d) 


S{t) = sinh ^ , C{t) = cosh ^ , 

T] Tj 


(23) 


with the hermitian matrix Xq = g(YqY^Y^‘^. The basic ansatz for the solution consists of fonr 
nnknown coefficients each for Q^^(t) and coefficients by implementing 

the Dirichlet bonndary conditions (0 1)Q^”^(1) = 0 and (1 O)Q^^(l) = 0 for the orbifold-odd 
prohles on the IR brane. Two more coefficients can be eliminated by imposing the Nenmann 
bonndary conditions {dt 0)Q^"^^(f)|t=i = 0 and (0 (9t)Q^^(t)|t=i = 0 for the orbifold-even 
prohles. Finally, the solutions for t E [1 — ? 7 ,1] are given by (Q = U, V) 










(24) 


where we use the short-hand notation 1 ,, = 1 — 77 . The fnnctions S{t),C{t) are given by (23) 
with Xg replaced by Xg = g{YYYqY/‘^. 

The derivation of the 5D qnark propagator in the region near the IR brane has been dis¬ 
cussed in detail in [29]. Appendix B contains all solutions that are relevant for the present work. 
Here we just comment that the basic solutions for the 5D propagator functions A^^(f,F; k"^) 
with A,B E {L, R} for f, t' E [1 — 77 ,1] are given in terms of the trigonometric fnnctions C{t) 
and S(t) in (23) with Xq replaced by 


s, = v/W+FIi, (25) 

where kE = kE/M^iK is the Enclidean momentnm normalized to the KK scale. The 77 depen¬ 
dence of the propagator enters only via the prodnct rjkE- As we will see below, this leads to 
a different behavior of the propagator depending on whether rjkE y^g or rjkE ^ y*g- 
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Calculation of the propagator functions 

It is instructive to discuss the calculation of the function in more detail, since it 

exhibits a sensitivity on the regulator r^, which is similar to that observed in the calculation of 
the loop-induced Higgs coupling to two gluons [29, 35-45]. Applying the (5-function regulator 
(19) and inserting the solutions for the external quark prohles (24) into (21), we obtain 




C{t) 

C{1,) 




Xd S{t) 




C(1 


g C(l^ 


'■RL 


{t, t'] k 


Q C(\^) 

Xd S{t' 


d,ll 

^RL 


{t, t'] k% 


m 

c(i,: 


g C{\ 


PrzVf(l^), (26) 


where C{t),S{t) are dehned in (23) with q = d. The propagator functions in the region near 
the IR brane for t, t' G [1,,, 1] can be found in Appendix B. 

We are not interested in the full dependence of V: since in the end of the 

calculation we will always remove the regulator (77 —)■ 0). However, since depends on 

the product rjkE via the 5D quark propagator functions and we integrate the function in (20) 
from zero to inhnite Euclidean momentum, we have to investigate whether the momentum 
integration commutes with the limit 77 —)■ 0. If we implement a momentum cutoff kE < Acut 
for the integral, the question can be reformulated as whether ( 20 ) yields the same results 
when imposing the constraints 77 3> 7/*n/Acut or 77 7/*n/Acut- Thus we need to investigate the 
ultra-violet (UV) behavior of T^^{ky) for large Euclidean momenta near the cutoff kE ~ Acut- 

Let us begin with the hrst scenario 77 3> y*v/Kcnt, where 77 is bounded from below. In 
fact, we also have to impose an upper bound 77 t/^u/Mkk, which is required in order to 
hnd reliable solutions for the 5D propagator functions in the region t,t' G [1,,, 1] [29]. When 
we consider large Euclidean momenta near the UV cutoff {kE ~ Acut), the allowed range 
of 77 implies the hierarchy kE 3> y^g/rj. Consequently, the function Sd = {X d + 
which is contained in the 5D propagator solutions, becomes approximately independent of the 
Yukawa-dependent term, such that Sd ^ ykE- In this limit, we hnd that {rjkE 3> y*g) 

TrlQ^'e) ~ (v^e) ^ (27) 


falls off with the third inverse power of the product rjkE- An analogous analysis for 
shows that it exhibits the same behavior as in (27). Since the imposed cutoff can be identihed 
with the effective UV cutoff of the theory near the IR brane, Acut ~ Atbv, the behavior in 
(27) refers to the case of a narrow bulk-Higgs scenario. 

We continue with the second scenario, where the (5-function regulator is bounded from 
above by 77 y^v/Acnt- This case represents the brane-localized Higgs scenario for Acut ~ 
Atcv- Consequently, the product rjkE is much smaller than y^g implying that Sd in (25) 
becomes approximately independent of the regulator, Sd ~ Xd- In this limit, we hnd {ykE 
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( 28 ) 


y*Q) 


_'n(2)t/i N f^ dtdt' W 2Xd Zd{k%) C{t') 

^Rd^E) Pl ^2 [^sinh2X, 1 + Z,(A;|)C(1,) 

, Xrf C(t') 


cothXd C(l^ 


- e{t - 1') X, 


S{t') \ C{t') 


/ 2X, Zdkl) S{t') 

\sinh 2 Xd 1 + Zd{k%) C(l^) 




C(l,j)/ C(l^) 


Pvivf(i,y 


cothXjC(i,) ' ' “c(i,);c(i,) 

Here we have introduced the structure [29] (for q = u,d) 

Z^{kl) = Y, R,CkE)YjRQCkE) , 

with the modihed Yukawa matrix Yg = (tanhXg/JYg)!^. We further need the ratio 


(29) 


RAik-E) 


- Ic^+i{ekE) I_CA + ^i^E) _ 
I-ca-I^^^e) Ica + \(^e) - IaA + \i.^kE) I_CA-lCkE) 


(30) 


of modihed Bessel functions, where cq, Cu^d are the bulk mass parameters of the 5D quark helds 
[4, 5]. In (28) we have not yet combined the terms inside the two round brackets. But when 
combining them we hnd that the t, f'-dependence completely cancels and the t, t' integrations 
become trivial, an analogous observation was made for the Higgs production process via gluon 
fusion in [29]. In addition the remaining ? 7 -dependence completely cancels. Finally, in the 
brane-localized Higgs scenario we hnd the results {rjkE y*Q) 




2Xd 


Zd{k 


2X. 


d 




[sinh 2 Xrf 1 + Zd{k%) 

2 Xd Zd{kl) 


sinh 2 Xd 

1 


+ 


-^d 


cosh^ Xd 
2X, 


(3)^ 


YdPuV^^’il 


(31) 




kj^ sinh2Xrf l + Zd{kl) Rq^e) sinh2X, 


d 


which are independent of the (^-function regulator. We have also included the hnal result for 
T^j^{k‘^), which can be obtained by an analogous calculation. For large Euclidean momenta 
kE 3> Mkk the structure Zq{k%) in (29) can be expanded as Zdi^k"^) ~ g'^YdYj + 0{k^'^). We 
observe that T'^dk^) reaches a non-zero plateau in this limit, which is in contrast with relation 
(27) valid in the narrow bulk-Higgs scenario. Consequently, the contribution of T)^i(/c|;) to 
the dipole coefficient (20) exhibits a dependence on the model under consideration. On the 
other hand, the function Tf^dk^) vanishes also in the narrow bulk-Higgs scenario and does 
not lead to a model-dependent contribution. 

Interestingly we could have obtained the same results for T^ik^ and T^Eiks) (31) if 
we had naively evaluated the extra-dimensional coordinates at t = F = 1“ instead of using 
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the regularized ^-function in (21). We have explicitly conhrmed that {r]kE y*g) 



^ 'T-)(2)t 

'V2 ^ 


V2 


(1 )-^d ^Rli^ 'i ^e)-^ d '^r\^ )) 

f'(1-) A?,„(l-, 1-; 4) M'g , 


(32) 


lead to the results (31). An analogous situation was encountered for the calculation of the 
propagator functions T±{k%) in the case of the Higgs production process via gluon fusion [29]. 


Final result for the Wilson coefficient 

The above analysis shows that the integrand of the dipole coefficient in (20) falls off with at 
least two inverse powers of Euclidean momenta ~ in the UV, which implies the finiteness 
of the integral. Thus, we are allowed to perform partial momentum integrations in (20) and 
find 


77, 8 g . 
^ 

AGpXt V 


lim IhMl 

ks^oo 4?7lfe 


1'°° dks kEm\ 

Jo {kl + mlf 


'^RLik's) ^E T^^{k%) \ 
mb k% + ml 2 Mkk /. 


(33) 


where all boundary terms at = 0 vanish. Based on the previous analysis only for large 
Euclidean momenta we can have a non-zero boundary term in case of the brane-localized 
Higgs scenario, where T^Likl;) approaches a non-zero plateau. This is accounted for by the 
first term in the outer bracket of (33). We can insert our results for T^Likl;) and T)j^(A;|;) in 
(31) into (33) and find 


77-89 




with 

9{X,,% 

9{X,,Y, 




2X. 


X 


4mb 

P 12 Zd{k 


sinh 2X. 


d Jo 


dkE kEml 
(/c| -f 


2',3 


2 X, 


mb 1 + Zd{k%) sinh2Xd 
, P+ kE Zd{k_ 


^ 2 ' 

''E, 


1 2X. 


(34) 


' 77 ( 3 ) 


2 k% + mll + Zd{kl) RQikE) sinh 2Xd 


vt’ii-) 


2X„ 




2Xr, 


brane 


__ 1 _ Y = +o^YY^Y + Oid^) 

Higgs sinh 2 X, 1 + g^YgY^ sinh 2 X, J, 




(35) 


narrow bulk-Higgs 


cosh"^ X, 


Y, = -q^Y,Y^Y, + 0{q^ 


The function g{Xq, Yq) is model dependent. To leading order in it only differs in the 

relative sign for a brane-localized and narrow bulk-Higgs. A similar observation was made 
for the KK tower contribution in case of Higgs production via gluon fusion [29]. We will see 
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numerically in Section (4.4) that this term emerges from the penguin diagrams exchanging 
KK quarks. Note that in this present paper we limit our analysis of the narrow bulk-Higgs 
model to the contributions involving the zero modes of the scalar doublet (see in particular 
Section 4.5). The contributions of scalar KK excitations have been studied in [20, 30]. 

We can generalize the results obtained in the brane-localized Higgs sector by allowing 
for two different Yukawa matrices and Y^'^ associated with orbifold-even and -odd quark 
prohles [41, 46]. In other words, we associate the correct-chirality Higgs coupling with Y^^ and 
the wrong-chirality coupling with Y^^, i.e. we replace in the Lagrangian for the Higgs coupling 
to down-type quarks Y^ —)■ Yj" in the hrst term and Yd —)■ Y/ in the second term of (18). We 
will refer to this model as the “type-H brane-Higgs” scenario [29]. We hnd that our previous 
analysis still holds, provided we use Y^ = (tank JCg/Xq)Y^^ for the modihed Yukawa matrix 
and Xg = We then obtain 




type-II 
brane Higgs 


2X, 2 Xd ^ 

sinh2Xd 1 g^YgYg^ sinh2Xd '' 




where to leading order in the KK contribution emerges from the correct-chirality 

Higgs coupling. At this order there is no difference between the original result (35) and (36). 
The wrong-chirality Higgs coupling only contributes at order 

In the following our paper concentrates on the RS model with a brane-localized Higgs 
sector and we set for simplicity Y^^ = Y^ = Yq. The only exception is Section 4.5, where we 
discuss some results in the narrow bulk-Higgs scenario. 


3.2 Gauge-boson contribution 

We continue with the diagrams (Ha) and (Hb) in Figure 1, where two internal quarks and one 
gauge boson are exchanged. The Wilson coefficients for the vector and scalar contributions 
are given by (with B = A,G,W, Z) 


K 


77.83 

B 


vector 


X 


27r / dkE I dtdt' BB{t'P b 

^LL (t-i t ) ^ 1 ) Pb Pl ^ ) 


16 ^ 16 ^ 48 


+ 



- Y, - 
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scalar 

77 , 8 s 


K, 


77,83 

B 


^ -2n dki 

4\/2Gi7’At Jo 


X 


U2 U3 UA 


32 


32 


dtdt' t; kl) ^t) Ve- (t) 


96 


A'^ 


l3 

TT 




A%L{t,t';k 


nib 


'^’VBt(t')vf(t') 


(37) 


where we introduced the matrices Pa = -Pg = l 2 x 2 and = (P+ + g'^/gf^P-), with gf = 
Tj — Qdsl; and gf^ = —Qdsl^. In case of the hh^-boson loop up-type 5D quark propagator 
functions (g = u) arise, otherwise we need to set q = d in (37). The quark propagator functions 
A^sit, t'k"^) are given explicitly in Appendix B for the brane-localized Higgs scenario. We 
remark that in case of the photon and gluon contributions to the Wilson coefficients {B = 
A, G) only KK resonances can contribute, therefore we have to subtract the zero mode 4D 
propagator (27rA;|;)“^ from Bsit',t', k"^) in (37). The structures^ ^±(f) and can be 

found in Appendix A. The coefficients are given by 


= ‘^Qd > 


K‘J = 2QdCFgi, n^ = Q. 


77 _^ 95 


— 2Qd [g'l) 


d\2 


2nr ’ 


gj/cj 

2nr 


80 ^2 


^83 _ ^ 
^ 27rr ’ 


4^ = 2(4)^ 


95 /' 


(38) 


27rr ’ 


where e is the 4D electromagnetic and the QCD 4D gauge coupling. The 5D gauge coupling 


of SU{2)l can be obtained from 


gj = 

27rr 


Amt 


— 4:\/2G 


1 _ 

ZiV^KK 


(1-4) + ^ 




which can be derived from the expansions of rfiw and v to leading order in given in 

the text below ( 6 ). Furthermore Qu = 2/3, Qd = —1/3 and Gp = {N^ — l)/{2Nc) = 4/3 with 
Nc = 3 being the color factor for quarks. The largest factors occur in case of the penguin 
diagrams exchanging KK gluons and hF^-boson modes. 

The scalar Wilson coefficient in (37) contains the propagator function F; —k"^/^), 

which is related to the 5D scalar propagator in general gauge via k) = 

_ i^«;aiar(^, _j^2 j j^g decompositiou (P = A, G, IF, Z), 


Rscalar., V. 2 n _ V- ^KK tt' ^tXn i^) OpXn ) 

^ ^ m\ ml — — it) 


(39) 


can be used to express it in terms of the vector-boson propagator ( 11 ) by means of the relation 


B 


scalar / 
B 


(t, t'] kp) — 


^KK JJJ_ 
e2 kl 


dtdf 


Bpit, t'] 0) — Bpit, t'] kp) 


(40) 


^The ± labels on the subscripts of Vg± (t) and Vg± (t) are only relevant for i? = lb and can be ignored 
otherwise. 
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We can use this equation to eliminate the brane-localized terms inside the structures Vg±(t) 
and Vg±(t) in case of the massive gauge bosons {B = W, Z). For example, the 5D Feynman 
rule for the DlUr vertex given in (A. 2 ) contains the term 



e 

t 



Lm^ 




(41) 


where = Y^P,2-Y}P21. The first term originates from the hfth component of the gauge- 
boson coupling to quarks, while the second brane-localized term is due to the Yukawa coupling 
of the IF^-Goldstone boson. We now insert (40) into (37) and perform partial integrations for 
the t, t' coordinates, taking into account that all terms on the boundary are orbifold-odd and 
therefore vanish. The partial integrations lead to derivatives acting on fermion prohles and 
propagators. We can use the equation of motions for the fermion prohles and the differential 
equations satished by the 5D propagators to show that all brane-localized terms contained 
in Vy^±{t) and V^y±(t) cancel. For example, the partial t-integration of the scalar Wilson 
coefficient in (37) leads to the term (for B = W) 


d. 




h2)t/ 


w- 


AU 

^LR 


{t, t'] k 


m. 




Ml 


KK 


Ml 


KK 


W 


AU 

^RR 


{t, t'] k 




RR{t, t ; /c|;) 


(42) 


where we have used that 


m^\t) 


t ; k 


Mkk 

T 7 ^ ) ^b) “ A4q(t) A|,^(t, t ; k\). 

tWKK 


(43) 


The last term in (42) cancels with the remaining contribution from the brane-localized term 
in (41). Furthermore, we will discard contributions that are suppressed by the strange-quark 
mass. 

In the last step we can perform partial integrations of the Euclidean momentum variable 
and hnally obtain [B = A,G,W, Z) 


fiB _ 
Sr, - 


K, 


77.83 

B 


A^/2GF\t 1 24 


^ ^fL-7^M + 27r rdkE [\tdt'vf^\t)PB 


X 


^lRit,t'-,kl) 


4 

(T/VN 


mb 


PBV^R^t') 


Qkp 

—^dk. 


3k% 2 

a 


(44) 


+ Al^(t,t';kl)PBVf\t') 



32 32 


k ^e) 


1 


where we have combined both the vector and scalar Wilson coefficients. The partial momentum 
integrations are required for numerical reasons, since momentum derivatives acting on fermion 
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propagators lead to complicated expressions that are very inefficient to evalnate. Note that 
the Wilson coefficients for B = A, G differ only in the factors see (38). Dne to 

the partial momentnm integrations we enconnter non-zero bonndary terms for large Enclidean 
momenta in (44), they are dehned by 


/?f^ = -27r lim kl 

ks^oo I 


= —27r lim fcl / dtdt'i?B(t', t; 0) ' (t) Pb 

ks^oo J 


( 2 ) 1 D 


(45) 




mb 


where q = u for B = W and g = d for P = A, G, Z. In case of the pengnin diagrams, in 
which photon (glnon) modes are exchanged, we have to subtract the zero-mode contribution 
from the full propagator function P^(t,t';0). The reason is that massless gauge bosons have 
constant prohles that lead to flavor-conserving interactions and therefore do not contribute to 
the Wilson coefficients. 


Calculation of the boundary terms R^l R^r 

In order to determine the boundary terms in (45) we need to know the UV behavior of the 
boson and fermion propagator functions, which is worked out in Appendix C. Using the results 
shown in equation (C.2), we can calculate the first boundary term in (45) and obtain 

Pf^ = 27r rdtPB(t,t;0)PfHt)P|Pf(t), (46) 


where we have to remember to subtract the zero-mode propagator in case of P = A, G, since 
only KK photons and KK gluons can contribute. We can further simplify (46) by using the 
explicit expressions for the propagator functions [10] 


BB(t, t] 0) 


BB{t, t] 0) 


1 

2 nfh% 


B = W,Z, 


1 



t^(l — 2Int) -I- 



B = A, G, 


(47) 


where B'^{t,t', 0) in the second line includes only the KK modes. Inserting (47) into (46) and 
applying the orthonormality condition for the fermion prohles (15), we obtain 


R 

R 

R 


A _ 
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W _ 
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Z _ 
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2 M^k 
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Dj23 
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(48) 


A _ {9^\ (^ 

V idiy) rhl 
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where = Rii- We recover the known overlap integrals 


)nn' 




fA' 


D Jnn' 




(49) 


{SD)nn'= I dtVP\t)P^V'£\t), ieD)nn'= I dtR Vf>\t) P.Vf > {t) , 


dn') 


Wt/ 


dn')i 


originally dehned in [10]. For the other boundary term in (45), we can use relation (C.9) and 
obtain 


RfR = 


27r 

Xh 


i?s(l-,l-;0)P®\l-)PB 


P_ q^Y}Y„ 


■q -*9 


1 + q^Y,Y^ 1 + q^Y^Y, 


12 


1 + q^y,y^ 


gYq - gY^ 


t 


21 


(3) /I -I 


1 + g^Y,Y, 


Pb 


(50) 




where Pij is a 2 x 2 matrix with zero entries except for the ij-component, which equals 1. We 
have omitted the terms at t = e, since the upper component of (t) and the lower component 
of {t) obey Dirichlet boundary conditions at the UV brane and therefore vanish. In order 
to obtain the last term we have used that the function vanishes for t' < t and 

we applied the equation of motion for the fermion prohles. We can further simplify (50) by 
performing a partial t-integration of the term involving dtSsit, t; 0 ) and by using the fermion 
equation of motions to show that 


Rm = 2 Rix 


jB 


2tt 


— Pb(i-,i-;0)p®Hi-)Pb 

Xb 


P^ 1 




P 1 




2 l + g^Y,Yg^ 


2 1 + g^Yg^Y, 


12 




21 


1 + " 1 + £)W,y, 


(51) 


f>Bcg’(i-). 


where we neglected a term suppressed by ms/rrib and where we recovered the term R^b- 
Applying the modihed boundary conditions of the quark profiles^, we obtain 
1 


pA _ pA 

^LR — 2 ^LL : 


dW _ ^ dW ^ ^ 'n(2)t/'i-i D ^ ^ ^uY^ ^ 

RlR - o RlL - -- Rl (1 ) -^12 ^P/J (1 ) , 
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TdZ _ ^ dX 
^LR — n^^LL 


2 'nf^ y/2mb 
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1 -f) Rm 


1 + 
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(52) 
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(3) /I -) 


1 + g^YdY„ 


x®!)'(i 


3ln the brane-localized Higgs scenario one can consistently calculate the fermion profiles by using modified 
boundary conditions at the IR brane without the notion of a regulator for the J-function [35]. Here we can 
apply {qY^ 1)p[P( 1-) = 0 and (1 - (1“) = 0. 
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where = -Rfjj. Yukawa-dependent terms appear in case of the massive gauge bosons and 
originate from the Goldstone degrees of freedom, which are localized at the IR brane. 

Finally, we have succeeded in obtaining expressions for the boundary terms (48) and (52), 
such that the complete Wilson coefficient in (44) can be evalnated nnmerically. The 

chirality-ffipped Wilson coefficients Cj.y^sg can be obtained from (44) and (45) by interchanging 
the label L ^ R. The bonndary terms can be calculated in analogy with the above steps, and 
we can use the results (46) and (51), for which we hnd = R^j^\l^r and R^j^ = Rfj^li-ir^R. 


3.3 Triple gauge-boson vertex contribution 

Finally we discuss the diagrams exchanging two internal gange bosons {B = W, G) and one 
qnark, see diagrams (Illa)-(IIId) and (IVa)-(IVd) in Figure 1. There are four diagrams each, 
involving vector and scalar components of the gange-boson propagators. We refrain from 
showing intermediate steps of the calcnlation bnt mention that we can proceed analogously 
as in the previous section and combine the vector- and scalar-boson contribntions. After some 
algebra we obtain the Wilson coefficients 
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32 
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96 




rUb 




dk 


^^E q2 


ks 


Bq^I , t] k^) \ , (53) 


where kww = 5'5/(2'7rr) and hgg = Note that the factor kqg for fhe triple glnon vertex 

diagram is larger by N"^ = 9 compared to k^q and comes with a relative sign. We recover the 
same bonndary terms that have already been calculated in Section 3.2 apart from constant 
factors. 


4 Analysis of the dipole coefficients 

4.1 Finiteness of the integrals 

In order to show the hniteness of the dipole coefficients in (20), (44) and (53) we need to know 
the UV behavior of the 5D boson and fermion propagators. We refer to Appendix C for the 
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corresponding derivations. For instance, the general behavior of the 5D (vector) gauge-boson 
propagator is given by (subscript B = A,G, W, Z) 




Vti’ 




(54) 


For large Euclidean momenta the propagator function is exponentially suppressed except for 
\t — t'\ ~ 1/kE- Integrating (54) along the coordinates t and t' we find 


"1 1 
dtdt'BB{t,t''-,k%) , (kE^l/e) 

kE 


(55) 


showing that the integral scales like k^'^ for Euclidean momenta kE 3> 1/e. Based on this 
analysis, and extending it to the case of the fermion propagator functions, we can formulate a 
power counting for integrals, where each extra-dimensional coordinate is integrated over the 
full interval. Excluding brane-localized terms, the counting in terms of Euclidean momenta 
can be formulated as 

A\b ^ {kE) Bw^z,a,g ^ {kE) [ dt ^ {kE) (56) 


where A, B E {L, R} and q = u,d for the quark propagator functions and with the additional 
condition that the last t-integration is not counted. This condition can be traced back to 
the conservation of the total 5-momentum. We can apply the power-counting scheme (56) 
to the penguin loops (44) and (53), showing that after the t,t' integrations the integrands 
fall off like k^^ for large Euclidean momenta. Thus the remaining momentum integration 
can be performed and yields a hnite result. This is in agreement with the Endings of [16], 
where the authors derived a power-counting scheme for the penguin diagrams treating the 
Yukawa interactions as small perturbations. The Higgs contribution contains two brane- 
localized vertices and our scheme (56) does not apply. In fact, the analysis of Section 3.1 shows 
that the propagator functions in the brane-localized Higgs scenario scale like T]^/j(/cl;) ~ k^^ 
andT|i(A;|) ~ const-|-(T(/cg^) for large Euclidean momenta. Since the Higgs boson propagator 
scales like k^'^ the Wilson coefficient is finite which is in agreement with the results of [16]. 

In summary, using our expressions for the 5D propagators with non-trivial boundary con¬ 
ditions at the IR brane we have confirmed the findings of [16] that the dipole coefficients are 
finite and calculable. This conclusion is also consistent with our results (57) derived in the 
KK-decomposed theory and discussed in the following section, where we can show that all 
dipole coefficients converge after summing up the KK towers. 


4.2 Connection with the KK-decomposed theory 

We can express the dipole coefficients, as defined via the amplitude (8), in terms of sums over 
zero-mode and KK-mode contributions. Starting from the expressions (20), (44) and (53) in 
the 5D framework we replace the 5D propagator functions by their corresponding KK repre¬ 
sentations. The appearing momentum integrals can be performed analytically and we obtain 
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the loop functions 13 , 4 (x) and h-ii{x), which are dehned via their integral representations 
(D.l) in Appendix D. We hnd the compact expressions 
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(57) 


where = ni^/ml, and g = n for i? = hh and g = d for i? = A, G, Z in the first line. The 
summation index m counts the contributions from the gauge-boson zero (m = 0 for the SM 
gauge-bosons) and KK modes (m > 1), while n counts the quark zero (n = 1, 2, 3 for the SM 
quarks) and KK modes (n = 4,9 for the first KK level and so on). We mention that there 
are no contributions from the massless zero modes (the SM photon and gluon), which implies 
that the summation starts with m = 1 in the first line for B = A, G and in the third line of 
(57). The ± superscripts on the overlap integrals and only relevant in the case 

of i? = W and can be ignored otherwise. The dehnitions of the overlap integrals can be found 
in Appendix A, while explicit expressions for the loop functions I 3 ,a{x) and Iq.ii{x) are given 
in (D.2). We note that when we insert the integral representations of the loop functions (D.l) 
into (57) we can identify the boundary terms and R^R) dehned in the 5D approach by 
(45), with the expressions {B = A, G, IT, Z) 
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in the KK-decomposed theory. Those terms originate from penguin diagrams where scalar 
components of the 5D gauge bosons are exchanged. In fact, we have also checked equation (57) 
by using the 4D Feynman rules listed in Appendix A and following the basic steps to obtain 
the dipole coefficients. The chirality-hipped coefficients Cy^sg can be obtained by replacing 

Kimfc ^ Vumk and {gi)nk ^ {gt)nk- 

We emphasize that there are two terms in each round bracket for the Wilson coefficients 
in (57). In the SM only diagrams with a chirality hip on the external 6 -quark line contribute 
to Cy, since the IT^ boson couples only to left-chiral quarks. Since in the RS model we can 
have also couplings to right-chiral quarks, there are additional contributions originating from 
diagrams where the chirality hip is performed on the internal quark line, which generates the 
factor rriq^/mb in front of the hrst term in each of the brackets in (57). When exchanging KK 
quarks in the loop this factor is large and enhances the contributions. 

We remark that we have numerically checked that all corrections in the RS model decouple 
with which is not directly obvious from the expressions (57). For instance, let us discuss 
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the KK contributions to C'^sg- stress that the loop functions can only take values 

from the compact intervals Iq{x) G [— 2, — |] and h{x) G [^,|], and are irrelevant for the 
discussion. We begin with the contribution of penguin diagrams that exchange SM quarks 
(n = 1,2,3) with KK W bosons {m > 1) in the first line of (57). Obviously the suppression 
by the squared KK hh-boson mass implies that the contribution decouples with Next, 

we discuss the contribution from exchanging the W boson (m = 0) with KK quarks (n > 
4). The overlap integrals scale like ~ ~ ~ ^kk n > 4, which implies 

that the second term in the round bracket of the hrst line in (57) decouples with M^. The 
hrst term in the round bracket is more subtle, since it is enhanced by the KK-quark mass 
niq^ ~ TTkk- However, numerically we observe that the summation over complete KK levels 
(n = 4,..., 9 for the hrst KK level and so on) leads to cancellations such that there appears an 
additional suppression. Hence, also the hrst term in the round bracket, when summed 
over complete KK levels, decouples with M^. In a similar fashion we can proceed with the 
contributions from the penguin diagrams with KK W bosons and KK quarks. The discussion 
can also be extended for the remaining Wilson coefficients. In fact the decoupling behavior 
with is apparent in the approximate expressions that will be given in Section 4.4. 

Finally, notice that the Wilson coefficients in the SM can be recovered from the second 
terms in and by summing only over the gauge-boson {m = 0 ) and quark zero 

modes (n = 1,2,3), and by replacing the overlap integrals with the CKM matrix elements 
Kr„s and with ^1,2,3 = «, c, t. 

4.3 Numerical evaluation 

The hrst step is to generate anarchic 5D Yukawa matrices, where each entry is bounded 
from above by y*, i.e. |(l(j)jj| < y^. The real and imaginary parts of each entry are random¬ 
ized with a hat distribution. If we expand the exact prohles of the W boson and the SM 
quarks in and only keep the leading terms, referred to as the zero-mode approxi¬ 

mation (ZMA) in [10], we can directly calculate the Wolfenstein parameters p and fj solely 
from the 5D Yukawa matrices. Next, we choose a random value for the bulk-mass parameter 
Ctjg G [—1/2,1], which corresponds to a localization of the right-chiral top quark near the IR 
brane. Working at leading order in we can determine the remaining eight bulk-mass 

parameters 2,3 from the experimental values for the six quark masses eval¬ 

uated at the scale p, = 1 TeV and from the Wolfenstein parameters A and A. Then, we choose 
a random value for Mkk G [1,10] TeV {Mg(i) G [2.45, 24.5] TeV) and calculate the whole set 
a^theo = {^u, nidi ^si nic, nib, nit, A, A, p, f)} using the exact expressions that are valid to all or¬ 
ders in Finally, we calculate the function x^{^) = X]n=i(^exp(^^) ~^theoin))"^/(T^,^p{n), 

where Xg^p contains the experimental values of the quark masses and Wolfenstein parameters 
with standard deviations given by ciexp- Points with x^(a;)/dof > 11.5/10, corresponding to 
less than 68 % CL, are rejected. Based on the procedure described above we generate six sets 
of 5000 RS points with different values for = 0.5, 1, 1.5, 2, 2.5 and 3. The upper value on 
I/* originates from requiring that the Yukawa sector remains in the perturbative regime [9]. 

We have implemented the integrals arising in the expressions for the dipole coefficients in 
(20), (44) and (53) in Mathematica. Since we expect the RS corrections to the SM Wilson 
coefficients to lie in the few percent range, we need to calculate the integrals to an accuracy of 
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Picture 

Average time to calculate 
Cj'y^sg for one RS point 

Average time fractions for each contribution 

^A/G /~iGG /~fh 

0^77,89 0/7.^ 0 ^ 77,89 0 ^ 77,89 0^89 '-^ 77,89 

5D 

571 min 

2% 4% 

21% 

32% 41% < 0.01% 

4D 

9 min 

16% 

46% 

38% < 0.1% 


Table 1: Time performance for calculating the Wilson coefficients and C^g in the 
5D and 4D (including 5 KK levels) pictures. The first column contains the average 
time needed in order to calculate the Wilson coefficients for one RS point on a 2.4 GHz 
Intel Core i5 processor, such that the results are compatible at the few per mille level 
in both approaches (see Figure 2 for more details). The additional columns show the 
relative fractions of time needed for the calculation of the six different contributions. 
Similar values are obtained in case of the chirality-flipped Wilson coefficients. 


a few per mille. We therefore set PrecisonGoal to 3 for the numerical integrations. Further¬ 
more we use a UV momentum cutoff such that kE < Acut = 100 Mkk, which improves the time 
performance without losing the required precision. It turns out that the numerical integra¬ 
tions over t and t' can be made faster by making the substitution t —)■ ^/tt = ln(t/e)/ ln(l/e) 
and analogously for F, which maps the integration region on the unit square. The first row of 
Table 1 compares the time performance of calculating the Wilson coefficients C'jry and C^g in 
the 5D and 4D pictures, averaged over many sets of RS parameter points. We need on average 
571 minutes per RS point on a 2.4 GHz Intel Gore i5 processor to calculate and C^g in 
the 5D approach. In more detail, the calculation splits into six parts belonging to different 
amplitude topologies. The least amount of time is required for the calculation of the Higgs 
contribution, since the t and t' integrations can be performed analytically, leaving over one 
momentum integral to be evaluated, see (20). Most of the computational time is needed for 
the KK contributions of the neutral gauge bosons, since the corresponding integrands involve 
all components of the 5D fermion propagator functions, in contrast to the hF^-boson Wilson 
coefficients in (44) and (53). 

We can compare our results from the 5D approach with the summation over zero and KK 
modes in the KK-decomposed theory, based on the results shown in (57). In order to achieve 
a consistency between both approaches at the few per mille level, we need to sum over five 
complete KK levels. The left plot in Figure 2 confirms that the results for calculated 
in the 5D and 4D approaches are consistent at the few per mille level. The results for the 
RS corrections relative to the SM Wilson coefficient are compatible at the 10% level in 
both pictures, as shown in the right plot in Figure 2. This presents a non-trivial cross-check 
of the formulas derived in Section 3. In both plots of Figure 2 we have focused on the real 
parts of but we have checked that the histograms look similar in case of the imaginary 
parts and also in case of C^g and the corresponding chirality-flipped Wilson coefficients. We 
need on average 9 minutes to calculate the KK quark and gauge-boson masses as well as the 
overlap integrals in the 4D formulation, see Table 1. Effectively there are only four different 
amplitude topologies, since and both depend on the same masses and overlap 

integrals, and analogously for C^gg and C^g^. Therefore we present combined time fractions 
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Re(C'77)5D 

Re(C7^)4D 


[%] 


Re[(C'7^)5D - 
Re[(C7^)4D - 


Figure 2: Compatibility of the results for the Wilson coefficient (left) and for 
the RS corrections relative to the SM Wilson coefficient ~ —0.20 (right) 

calculated in the 5D and 4D (including 5 KK levels) pictures. Both histograms contain 
RS parameter points with different values for and Mkk £ [1,10] TeV. The vertical 
dashed lines denote the median values of the corresponding distributions. 


for those Wilson coefficients in Table 1. When we require a consistency of a few per mille for 
the calculation of the Wilson coefficients between the 4D and 5D approaches, we find that the 
summation over KK levels is faster by a factor of order 60. As a consequence, after we have 
verihed that the results in the 5D and 4D approaches agree at the required level of precision, 
we will implement the equations in (57) for the numerical calculation of the Wilson coefficients 
for most of the RS points used in the phenomenological analysis in Section 5. 

4.4 Approximate expressions 

We emphasize that the integral expressions (20), (44) and (53) and the corresponding results 
in the KK-decomposed theory (57) are formally valid to all orders in All numerical 

results that will be presented are calculated from those equations according to the procedure 
described in (4.3). However, in order to better understand the size of the different contribu¬ 
tions from the diagrams in Figure 1 we will also derive some approximate formulas in this 
section. The hrst step is to parametrize the RS corrections relative to the SM Wilson coeffi¬ 
cients by 

C^^’°(p.vv) = , 

+ C^7T''(Rkk) + ^^7''/''(Rkk) (59) 

+ ^77 IhKKj + <-^77 IhKKj , 

where we distinguish corrections that arise from the exchange of only zero modes 

and of loops including at least one virtual KK particle C^®’^^(p.kk)- The individual zero-mode 
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RS corrections 

Median values of the distributions in [%] 
lU WW Z A G GG h 

\c^\f^w)\/EB\c^'M\ 

\C!g^^^w)\/EB\Csg\^^w)\ 

34.8 65.0 0.07 - - - 0.001 

99.6 - 0.4 - - - 0.007 

|C^7T''(/^Kk)|/Eb|C^7T''(/^Kk)| 

33.9 51.2 7.7 0.0001 0.1 - 7.0 

52.8 - 24.1 0.0004 0.005 0.6 22.0 


Table 2: Median values of the distributions in the left column based on RS points with 
= 3 and Mkk G [1,10] TeV. The median values can be used to estimate the relative 
size of the RS corrections arising from the exchange of only zero modes 
and of loops including at least one virtual KK particle Similar values are 

obtained in case of the chirality-flipped Wilson coefficients. 


contributions for B = W, WW, Z, h are dehned at the electroweak scale yw ~ rnw 

and are given simply by setting m = 0 and summing over n = 1,2,3 in (57). On the other 
hand the KK contributions 0^’^^(/iKK) are defined at the KK scale pkk ~ IbfKK- Analo¬ 
gous parametrizations hold for the chirality flipped Wilson coefficients and for C^g’^{yw) and 
^rs,kK(/^^k)^ where the triple gluon-vertex contribution must be included. In the SM the 
contribution to the Wilson coefficients at leading order is given by the penguin diagrams (II) 
and (III) in Figure 1, in which virtual bosons and up-type quarks are exchanged. The 
charm- and top-penguin diagrams yield comparable contributions, since the product of the 
CKM matrix elements are of similar size, |At| ~ |Ac|, where \q = V*gVgb. Making use of the 
unitarity of the CKM matrix, which implies A,j -|- Ac -l- A* = 0, the SM Wilson coefficients at 
the electroweak scale are given by 

-1) + h{xy) - T, ^ - 1, (60) 

where Xt = The loop functions h^x) and Igi^x) can be found in Appendix D. Using 

^^t,poie = 176.7l3;4GeV [47] we hnd Cy^{yw) ~ —0.20 and Cf^(pvv) ~ —0.098. 

In a first step, we look at the relative size of the RS corrections and (/^kk) 

based on a set of RS parameter points. To this end, we compare the median values of the 
distributions obtained from calculating |C'^’g^(/ipi/)| and |C'^’™(/iKK)| and normalizing them 
to the total sum of each (absolute) correction \^ 7 'y^ 8 gif^w)\ and I^^’,8 ^(/^kk)|- Table 2 
shows the results. The general pattern is that the penguin loop diagrams with hF^-boson 
exchange give the largest corrections, which is also true for different values of y^. In case of the 
zero-mode contribution gg we find that the largest corrections are given by the deviations of 
the overlap integrals and in the RS model with respect to the CKM matrix elements 
F)* and Vtb in the SM, and by the coupling of the W boson to right-chiral quarks. Those 
corrections stem from the non-flatness of the lU-boson prohle and from deviations of the 
exact (Z 2 -even) quark prohles from the ZMA expressions [10]. The zero-mode contributions 
from the Z and Higgs bosons arise due to their flavor-changing couplings to quarks in the RS 
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model, but they are suppressed by small down-type quark masses and can be neglected. To 
leading order in /M^-^ we find 
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(61) 


with the function F{c) = sgn(cos ttc) a/(1 -|- 2c)/(I — e^+^'^) [10]. The three-component vectors 
and for n = 1, 2, 3 form the columns of the unitary matrices L7 and which define 
the CKM matrix Vckm = Ud in the ZMA [10]. The corrections in the first line of (61) stem 
from the W^undR coupling, while the remaining terms include corrections to the CKM matrix 
elements Vtb and in the RS model. Concerning the KK contributions we find, contrary to 
the observation made in [17], (independently of y^) that the triple gluon vertex contribution 
is subdominant and does not enhance the chromomagnetic dipole coefficients. In general we 
find that the penguin diagrams with the exchange of photon and gluon KK modes yield very 
small corrections. At last we can compare the relative magnitude between and 

Numerically, we find that both contributions are similar in size for ^ 2. For 
larger values of y^ the KK contributions dominate in size. 

In the next step, we will take a closer look at the main KK contributions to the Wilson 
coefficients and derive approximate expressions. 


W- and Z^-boson contributions 


We begin with the KK contribution of the W^- and Z-boson penguin diagrams to the dipole 
coefficients. The dominant contributions come from diagrams in which charged/neutral scalar 
zero modes (stemming from the fifth component of the 5D gauge-boson field and the Goldstone 
bosons in the Higgs sector) and KK quarks are exchanged, which are implicitly included in 
the first two expressions in (57) for the Wilson coefficients.^ In this case we are allowed to 
take the limits 3> 1 and 1 for the loop functions in (57), leading to Iq{x) ~ —1/2, 

hix) K. 5/12, Is{x) ^ —1 and ~ 5/12. We find 
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^If we would not include the contributions from the Goldstone bosons, the diagrams with gauge-boson 
zero-modes and KK quarks would be suppressed to leading order by In this case the contributions 

from KK gauge-bosons and KK quarks would be dominant, since they contribute already at order 
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where the boundary terms and Rfn are given in (48) and (52). Since the limits of the loop 
functions we have taken are not valid in case of SM quarks, we have to subtract the contribu¬ 
tions from the quark zero modes in (62). We observe that the corrections of C^’^^(/iKK) and 
^ww,kK(/^kk) add up constructively, since = Qu gfH^Tir) and = gl/{2TTr). We can 
further simplify the boundary terms RYi^ and RYr in (62) by expanding them in 
and neglecting terms that are suppressed by mg/nib. We obtain 
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where in the case of RYr and Rrr we have implemented the relation [35] 
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We checked numerically for different values of that the approximate formulas (62) together 
with (63) are accurate at the 10% level compared with the exact expressions. We emphasize 
that the approximate expressions are independent of the masses and prohles of the KK quark 
and gauge-boson modes. In (63) we encounter terms including products of three Yukawa 
matrices UVJU and YdY^Y^ originating from the IR brane-localized terms in and R^r 
in (52). They originate from diagrams exchanging and Z Goldstone bosons with a chirality 
flip on the internal KK quark line as shown in Figure 3. Those terms yield the dominant KK 
corrections for not too small values of the Yukawa matrix entries, which is approximately 
fulhlled for RS points with > 1. In this case we can derive simpler expressions and hnd 
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Figure 3: For y* > 1 those diagrams give the main KK corrections for 

the transitions b ^ s'y and b ^ sg at the one-loop level. Internal solid lines labelled 
by Ur denotes the exchange of singlet up-type KK quarks, while Ul,Dl imply the 
exchange of SU{2)l doublet KK quarks. Crosses denote a chirality flip on the internal 
KK quark lines. Here, dashed lines labelled with W or Z denote the contributions 
from the corresponding Goldstone bosons in the Higgs sector. 


where the chromomagnetic dipole Wilson coefficients can be obtained from the 

hrst (second) line in (65) by sending Qu —)■ 1 (Qd —^ !)• Moreover, gjyg^ analogously 

by the expression in the hrst line of (65) with set to 1. We checked numerically that the 
approximate expressions are valid at the 10% level with respect to the exact expressions for 
RS points with y* > 1. Approximate formulas for the chirality hipped Wilson coefficients 
^^’™(hKK) can be obtained from (65) by making the replacements L -H- i?, -H- , 

T/W- V yw- yw+ V T / n /+ 

^203 ^ '"^203 ^303 —?■ ^393 . 


Higgs contribution 

The diagrams contributing to (^^.(^(/iKK) involve the exchange of the Higgs boson with KK 
quark modes. For the exchange of KK quarks we can use that x'l" 3> 1, allowing us to take the 
limits hix) = t/{2x) + 0(x~^) and l 4 (x) = l/(12a;) -|- 0(x~^). The contribution associated 
with Islx) dominates, since this loop function is less suppressed than l 4 (x) in the considered 
limit and the contribution is enhanced by md„/mb, which is a large factor for KK modes. If 
we only keep the corresponding contribntion associated with Islx), we obtain approximately 
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( 66 ) 


where we have expanded the expression to leading order in and neglected ms/rrib- 

snppressed terms. The corresponding expression for is given by (66) with Qd —)■ 1. Nn- 

merically we have checked that (66) is accnrate at the 10% level with respect to the exact 
expressions. Note that the YdY^Yd structure in (66) originates from the leading order expan¬ 
sion in of the function g{Xd, Yd) dehned in (35) and gives the dominant contribntion 

for > 1. In fact, this term exactly cancels the expression C^^^(/iKK) in (65). Conseqnently, 
fo^^ y* ^ 1 fhs KK corrections from the Z-Goldstone boson and Higgs diagrams cancel to very 
good approximation. 
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Dependence on M^^i) and 

Figure 4 shows histograms of the (absolute) KK corrections |C'^® 8 ™(/iKK)| and |(F^^ 8 ™(hKK)| 
for a set of RS parameter points with y* = 3 and M^(i) = 10 TeV. We choose y* = 3 to 
obtain maximal effects, while still staying in the perturbative regime for the Yukawa sector. 
The value Mg(i) = 10 TeV is close to the lowest KK gluon mass that is consistent with the 
tree-level analysis of electroweak precision data. The distributions can be well described by 
the approximate formulas given in (65) and ( 66 ). For different values of M^(i) and ?/*, the 
corresponding distributions can be obtained by the formula 
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which is a good approximation for y* ^ 1. For smaller values of the KK contributions do 
not follow a simple scaling law with ?/*. An analogous equation holds for the distributions of 
the chirality-flipped Wilson coefficients. In order to get a rough estimate for the typical size of 
the KK corrections, we can calculate the median values of the distributions of |C'^^8s^(hKK)| 
and hnd (y* > 1) 


Median (^|C'f^®^™(/iKK)|) ~ 077,89 X x , ( 68 ) 

with = 0.012 and asg = 0.0073. In case of the median values of |(F^® 8 ™(/iKK) | the coeffi¬ 
cients read 0 , 7 ^ = 0.020 and dsg = 0.012. These coefficients represent the median values of the 
distributions shown in Figure 4. We observe that the KK corrections to the chromomagnetic 
dipole coefficients are (approximately) smaller by the factor k^/{k^ + ~ with 

respect to the electromagnetic dipole coefficients. Furthermore, we hnd the general pattern 
that the chirality-hipped Wilson coefficients are enhanced, which can be explained by the 
diherent localization of the left- and right-handed bottom-quark prohles. The left-handed 
bottom quark prohle, which enters and (^^^’^^(/iKK), is more localized towards 

the IR brane (cf,^ = cq^ > Cb^ = Cd^) and is thus more sensitive to havor-violating ehects. 
This hierarchy of the bulk mass parameters is due to the large mass diherence of the top and 
the bottom qnark, which reqnires that F{cb^) > F{cbj^). 

We can (approximately) relate our results with the numerical analysis of the Wilson coef- 
hcients performed in [17], where the case of = 3 and Mg(i) = 2.5 TeV was discussed. When 
we consider M^(i) = 2.5 TeV we hnd that the corrections C^®’^^(Mkk) are larger by a factor 
of ronghly 5 compared with [17]. In case of C^®’^^(Mkk) we hnd that the corrections are 
similar in size. Concerning the corrections to the chirality-hipped Wilson coefficients, we hnd 
that they are larger by a factor of ~ 2 with respect to Cj^^sg, while [17] reported a stronger 
enhancement by one order of magnitude. While we have been unable to trace the origin of 
these discrepancies, the fact that we have performed our analysis using both the 5D and 4D 
formulations of the RS model and fonnd consistent results in both approaches provides a 
highly non-trivial cross-check of our calculations. 
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\^8g (A'-Kk)! \^8g (/^Kk)| 


Figure 4: Absolute corrections from KK modes to the Wilson coefficients at the KK 
scale for a set of RS points with = 3 and M^(i) = 10 TeV. The vertical dashed lines 
denote the median values of the corresponding distributions. The size of the corrections 
are equally distributed among the real and imaginary parts of the Wilson coefficients. 


4.5 Comment on the narrow bulk-Higgs scenario 

We have observed that the sum of the KK contributions and cancels to 

a very good approximation for y* ^ 1. In this section we investigate whether this cancellation 
still holds in the narrow bulk-Higgs model. 

The Higgs contribution was already calculated in Section 3.1, and the hnal result has been 
given in (34) and (35), including the case of a narrow bulk-Higgs. It follows that the first term 
(containing the YdY^Yd structure) in the approximative formula for (pkk) in (66) must 

be multiplied with a minus sign in the case of a narrow bulk-Higgs. Concerning the Z-boson 
contribution we focus on the scalar diagram (Hb). For ?/* > 1 the dominant corrections are 
due to the exchange of the Goldstone Z boson. The corresponding terms can be extracted 
from in (37). We can proceed analogously to the calculation of the Higgs contribution 

discussed in Section 3.1. Therefore, we refrain from giving more details here and quote the 
hnal result (?/* > 1) 
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The scalar Z-boson propagator behaves like 1~; k%) = L^m|/(27re^/c|;MKK)+^(^E^) 

for large Euclidean momenta, rendering the integral hnite. The function h{Xq, Yq) in the hrst 
line of (69) is given by 
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U = jY,Y,% + 0(e' 


where the difference between the brane-localized and narrow bulk-Higgs scenario is to leading 
order the relative factor —1/3. Thus, in the narrow bulk-Higgs scenario the approximate 
expression for (^^^(pkk) in (65) must be multiplied with —1/3. Finally, adding (^^/^(Pkk) 
and C'^™(pkk) we obtain approximately (for ?/* > 1) 
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The corresponding expression for the coefficient of the chromomagnetic dipole operator is 
obtained by replacing Qd —)■ 1. The structure YdY^Yd cancels in the brane-localized Higgs case, 
while there remains a non-zero contribution in case of the narrow bulk-Higgs scenario. This 
observation and the factors 0 and in (71) were hrst encountered in [48, 50] for the case 
of lepton penguin loops. In fact, we can exactly reproduce the result (26) in [50] for the 
Higgs contribution in the lepton sector from equation (71) by replacing Qd ^ Qe = 
and by accounting for factors in the dehnition of the Wilson coefficient. We note that while 
the contributions from the neutral scalars cancel for y* > 1 in the brane-localized Higgs 
scenario, we still have left over the (dominant) contributions from the charged Goldstone 
bosons. The latter contribution is absent in case of the leptonic dipole coefficient for the 
transition k —)■ Ij'-j in the minimal RS model, which does not include right-chiral SU{2)l 
singlet neutrinos. However, a non-zero contribution from neutral scalars would be present in 
case of the RS model with custodial protection, which can be found in [20, 49]. 

Finally, we remark that in order to calculate the contribution of the charged Goldstone- 

bosons in the narrow bulk-Higgs scenario, we need to perform t, t' integrations over matrix¬ 
valued functions mixing Yu with Yd- Since we could not handle those integrations in a semi- 
analytic way we will therefore conhne our analysis to the brane-localized Higgs scenario in the 
remainder of this paper. 
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4.6 Renormalization-group running to the meson scale 

In the previous section we have analyzed the corrections to the SM Wilson coefficients from the 
zero modes dehned at the electroweak scale nw ~ raw and from the KK particles 

at the KK scale /xkk ~ d^KK- For the phenomenology we are interested in the 
Wilson coefficients Cj^^sgifJ'b) at the meson scale j^b ~ ^b- When running down from higher 
scales down to /Xb, QCD effects generically lead to a mixing between dimension-6 operators. 
The general effective Lagrangian for a new-physics model at a high scale (/iKK in our case) 
can be written in the form 


r - 9 ^\ 
CeS ^ A, 


E E 

L q=u,c,t 2=1,2 


E E E ^iq),LA ^ ^id),LA 

2=1,2 A=L,R ^ q=u,c^t,d,s^b 


^ Q, C ^ C, L ^ R 


(72) 


where we adopt the notation of [51]. Here are the charged current-current operators, Qr^y^sg 
are the dipole operators, and (5i,2(^,-B), Qi, 2 A,B) are neutral current-current operators 
(including the four-quark QCD and electroweak penguin operators of the SM). In the RS 
model such operators are induced by the exchange of the heavy KK modes of the Z boson, 
photon and gluon. For simplicity, however, they will be neglected in our analysis. When 
running down from /ixK ~ Mkk to fiw ~ we consider only the mixing between Qjy 
and Qsg, which accounts for the dominant evolution effects. Between the electroweak scale 
and the meson scale we further include the RS contribution to the charged current-current 
operator = 4(sj7^PLCj) (ci'y^PLbi), which is also important in the SM calculation. The 
corresponding Wilson coefficient in the SM reads C 2 ^^^{nw) = —A/^t ~ 1- In the RS model 
the hF-boson coupling to quarks receives corrections, which we include later in the analysis 
by dehning the Wilson coefficient 
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The overlap integrals dehned in relation (A.6) of Appendix A. 

Let us outline the basic steps how we evolve the Wilson coefficients down to the meson scale. 
We need the evolution matrix U{fii, ^ 2 ) which can be calculated from the anomalous-dimension 
matrix 7 of our operator basis, which is contained in the more general basis considered in [51]. 
The running between scales is accomplished at leading order by 
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Figure 5: Corrections to the effective Wilson coefficients in the SM ~ —0.30 

and C|^ ~ 0.15 at the 5-meson scale /i^ = 4.8 GeV for a set of RS points with y* = 3 
and Mg(i) = 10 TeV. The vertical dashed lines denote the lower (25%) and upper (75%) 
quartiles of the distributions. The distributions look similar for the imaginary parts of 
the Wilson coefficients. 


where 7^°^ includes the eigenvectors of the transposed anomalous-dimension matrix The 

matrices V diagonalize 7^°^ , such that V is diagonal. Note that 7*^°^ and /do = 

(33 — 2/)/3 depend on the number of active flavors /. Between the scales //kk and Hb we 
integrate out the top quark, such that the evolution matrix splits into two parts, 

u (/ib, /iKx) = (p-b, Pvv) (pvV; h-Kx) • (75) 

The RS corrections at the KK scale, coming from integrating out heavy KK resonances, are 
contained in the coefficient C'^®’^^(pkk)- The evolution down to the electroweak scale is 
given by Fkk) <?^^’^^(/^kk)- At the electroweak scale the W boson 

and the top quark are integrated out. Matching on this new effective Lagrangian we include 
the contributions from the boson and fermion zero modes, which are given by = 

(T®^(pvv) + where C^®’°(pvv) contains the zero-mode corrections to the SM co¬ 

efficient. Next we evolve this contribution down to the meson scale. The effective Wilson 
coefficient reads 

C(pb) = fd®^(pb) + U{fXb, I^kk) + U{fib, fj^w) , (76) 
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Figure 6: Left (right) is shown the approximate linear correlation between the RS 
corrections to the (flipped) electro- and chromomagnetic dipole coefficients for RS 
points with = 3(0.5) and Mg{i) = lOTeV. A similar correlation is found for the 
imaginary parts of the Wilson coefficients. There is no correlation between C^^gg{fib) 
and their chirality-flipped counterparts. 


where the SM Wilson coefficients are given by ~ —0.30 and Cg^{fib) ~ —0.15. 

Performing all steps including the dipole and the charged current-current operators, the RS 
corrections to the electro- and chromomagnetic dipole operators at the R-meson scale are 
given by 

= 0.475 C^®’™(pkk) + 0.123 C^®’™(/iKK) + 0.667 
+ 0.092 C'8''/’°(/iiy) - 0.174 , (77) 

Cf-^ifib) = 0.522 + 0.702 - 0.080 • 

The numbers in front of the KK corrections have been calculated for pkk = 

1 TeV, Hb = 4.8 GeV and jj^w = 80.4 GeV. In our numerical analysis we set pkk = ATkk 
for each RS point. Relation (77) also holds for the chirality-flipped Wilson coefficients, since 
(massless) QGD is blind to the fermion chirality. 

Figure 5 shows the RS corrections to the SM values of the dipole coefficients at the B- 
meson scale for RS points with y* = 3 and Mg(i) = 10 TeV. In general, the RS corrections to 
the SM Wilson coefficients Cj^{y.b) and G|j^(pb) lie in the few percent region. On the other 
hand, the relative corrections to the chirality-flipped Wilson coefficients are large, since in the 
SM G7.^(/ib) and Gf^(/ib) are suppressed by mg/mb. The dominant contributions to C^^{y.b) 
are given by the RG-evolved KK and zero-mode corrections G^®’™(pkk) and G^®’°(/ivi/)- An 
analogous statement holds for G|^®(pb). The mixing of the chromomagnetic dipole Wilson 
coefficients into Cf^{fib) yields a correction of roughly 10% for y* = 3 and 7% for = 0.5. 
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Finally, we remark that C^^{nb) and are linearly correlated, which can be seen in 

Figure 6. This is expected, since the main contributions arise from penguin diagrams contain¬ 
ing VF^-boson modes, and they only differ by the factor = l + Qu = 5/3. The 

coefficients and as well as and Cfg{fib), are however largely uncor¬ 

related. 


5 Phenomenology 

5.1 Branching ratio Br(.B -X^sT) 

We begin with the CP- and isospin-averaged B —)■ branching ratio, which is one of the 

cleanest observables in B physics from a theoretical point of view. Measurements lead to the 
combined result Br(i? —)■ Xs7)exp = (3.43 ± 0.21 ± 0.07) x 10“^ [54] for the branching ratio 
dehne with a lower cut > Eq = 1.6 GeV on the photon energy in the meson rest frame. The 
SM prediction at NNLO reads Br(ii —)■ Xs 7 )sm = (3.36 ± 0.23) x 10“^ [55] for Eq = 1.6 GeV, 
showing that both values are compatible at the la level. In order to estimate the effects of 
the RS model we use the approximate formula (for E^ > 1.6 GeV) [56] 

Bt{B ^ X,7) = Br(R ^ X,7 )sm + 0.00247 - 0.706, (78) 

where the Wilson coefficients have to be evaluated at the R-meson scale fib = 4.8 GeV. While 
all known non-perturbative contributions (see in particular [57] for an estimate of non-local 
hadronic effects) are taken into account, the RS corrections are included at leading order in 
tts- This is accurate enough to estimate their impact, because these effects are generally small. 

The dominant corrections in (78) stem from the last term in the round bracket, which is 
proportional to ReG^®(/if,). The squared contributions (and in particular the chirality-flipped 
Wilson coefficient) have only a minor impact. Since the KK contributions are approximately 
proportional to y^, the biggest effects can be expected for large values of y^. There exists an 
upper limit y* < y^ax when requiring that the Yukawa sector remains in the perturbative 
regime, and it is conventional to choose the value y^ax ~ 3 [9]. We have generated RS param¬ 
eter sets for different values of y^ and Mg(i) (Yukawa matrices and quark bulk masses), which 
correctly reproduce the SM quark masses and the Wolfenstein parameters, see Section 4.3 
for more details. In the left plot in Figure 7, we show predictions for the branching ratios 
Br(R—)■ Xs'-f) and Br(R —)■ which will be discussed in the next section, for a large 

set of RS model points with Mg(i) = 10 TeV. The black (green) points are obtained with 
y* = 3 (0.5). We find that more than 90% (99%) of these points he within the experimental 2a 
bands. The RS corrections to Br(R —)■ X^y) approximately scale with in the region where 
y* > 2. For smaller values of ?/* there is no simple scaling dependence. In general, we hnd that 
the size of the RS corrections to Br(R —)■ X^y) is strongly dependent on y^,, in contrast to the 
observation of [17], where no signihcant correlation in their numerical scan was reported. If 
we require that at least 10% of the RS points he within the 2a error margin, we can derive 
the lower bound Mg(i) > 3.4 TeV for y^ = 3. This bound cannot compete with the constraints 
from a tree-level analysis of electroweak precision data at 95% GL. On the other hand, if we 
set Mg{i) = 2.5 TeV, which is the lowest value allowed from the direct search of resonances in 
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Br{B ^ X, 7 ) [10-4] Bv{B ^ Xsj) [IQ-^j 

Figure 7: Left is shown the branching ratio Br(i? —)■ Xgl~^l~) with respect to the 
inclusive radiative decay Br(i? —)■ Xg'j). The right plot shows the time-dependent CP 
asymmetry Sk*-/ as a function of the branching ratio for the decay B —)■ In both 

plots the light gray and blue bands show the la experimental error margins while the 
area between the dashed lines contains the SM prediction with Icr uncertainty. All black 
(green) points represent possible RS scenarios with = 3 (0.5) and Mg(i) = 10 TeV. 


the invariant mass spectrum of tt production by the ATLAS [58] and CMS [59] collaborations, 
the maximal Yukawa value can be constrained from above to y* <2. 

Let us comment on two further constraints on the RS parameter space. First, we consider 
the CP-violating observable ex in kaon mixing, which can receive large corrections in the RS 
model due to a strong chiral enhancement of the four-quark operator (^4 = {dRSL){dLSR), 
after performing the RG running from Mkk down to the kaon mass. When we impose the 
constraint that the RS prediction for ex lies in the 2 ct region of the SM prediction we find that 
roughly 15% (0.7%) of the black (green) points in Figure 7 survive. The fraction of allowed 
points decreases with smaller values of y*, since the RS corrections to ex are approximately 
proportional to l/y"^. Still, the shape of the distribution of points is not strongly affected, 
since ex is uncorrelated with the observables discussed in this paper. Secondly, we can discuss 
the impact of Higgs physics, where the strongest bounds arise from the signal rates of the 
Higgs decaying into pairs of electroweak gauge bosons. Comparing with LHC data one finds 
the condition Mg(i) > (15 — 20) TeV x (|/*/3) at 95% CL [27]. Applying this bound to the RS 
points with Mg(i) = 10 TeV would exclude the black points (?/* = 3) but still allow for the 
green points (?/* = 0.5). 

5.2 Branching ratio Br(.B — Xsl^l~) 

Next we consider the inclusive decay B ^ Xg l^l~ in the low region of the dilepton invariant 
mass, 1 GeV^ < < 6GeV^. From the latest measurements of the branching ratio in the 

low dilepton mass region from Belle Br(i? —)■ Xg = (14.93 ± h.Odlg;^!) x 10“^ [60] 

and Barbar Bt{B —)■ = (IG.Olg;^ ± 1.8) x 10“’^ [61] we take the combined 

value Br(R —)■ Xsl~^l~)eyr.p = (15.8±3.7) x 10“’^ [62], which is in good agreement with the SM 
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prediction Br(5 —)■ Xsl~^l~)gj^ = (16.2 ±0.9) x 10“^ [62]. For the calculation of the branching 
ratio in the RS model we need to take into account the electroweak (leptonic) penguin operators 
Qgpo and Qgpo by adding the effective Lagrangian Ces = ^*=9,10 (CiQi + CiQi) to (72). 

The operators are dehned by 

2 2 
<29 = ^ {s7i,PLb){lYl ), {s7i,PLb){h^75l ), (79) 

where as always the chirality-flipped operators can be obtained by replacing —)■ Pr- In the 
SM Cgi^fiw) and Cio{fiw) are loop suppressed. In the RS model corrections are induced at 
tree-level due to the flavor-changing couplings of the Z boson, the Higgs and the Z-boson and 
photon KK modes. The Higgs contributions are suppressed by lepton masses and can therefore 
be neglected. To leading order in the RS corrections to the Wilson coefficients are 

given by [13] 



where Oe = Qd = -1/3, Qi = -1, 5^/ = -1/2 ±s^ and = si,. The functions (Ad) 23, 

(A'£))23 and (5/))23 are dehned in (49), while (Ad)23 and (A'^)23 can be found in [10]. Figure 8 
shows the distributions of these Wilson coefficients at the R-meson scale. The coefficients 
are much larger than their chirality-hipped counterparts since the left-handed b- 
quark is more localized towards the IR brane than the right-handed one, F{cb^) 3> F{cbj^). 

We emphasize that the new-physics contribution (/i^) is largest for RS parameter points 
with small values of y^. The reason is that the overlap integral (A/))23 in (80) can be expressed 
approximately by {Ad )23 ~ F{cs^) F{cb^) [10], showing its sensitivity on the localization of 
the left-handed strange- and bottom-quark prohles. For smaller values of y*, the prohles are 
shifted towards the IR brane, in order to reproduce the correct bottom and strange quark 
masses, and the overlap integral increases in magnitude. For the branching ratio in the low 
G [1, 6] GeV^ region, we have implemented formula (3.9) in [63] and hnd 

Br(R ^ X, /+/-) = Br(R ^ X, /+/")sm 

± 10"^ X [l.41 Re - 0.74 Re ± 2.81 Re - 0.059 Re - 4.65 Re Gfo® 

± 0.074 ReGfo® ± 30.18 (1^7^^^ + l4^7T) + 0-52 (81) 

± 0.52 (|GfoY + \Cio\^) + - 0.008Re(G9^^G9^^* ± Gfo^Gfo®*) 

± 2.42Re(G7^^®G9^^* ± 4^®*) - 0.017Re(G7^^®G9^®* ± , 
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where all Wilson coefficients have to be evaluated at the 5-nieson scale f^h■ In our analysis 
we have used C7^(/ife) —0.30, 4.07 and ~ —4.31. We note that the 

coefficient of Re is smaller than naively expected. The reason is that the corresponding 

coefficient results from two terms that are interfering destructively, and the difference is rather 
sensitive to the SM values of the Wilson coefficients and Cg'^(/ife). We stress, however, 

that this sensitivity does not have a large impact on our analysis, since the corrections to the 
branching ratio in (81) are dominated by the RS corrections to Cio(/ife). 

The left plot in Figure 7 shows the branching ratio Br(.B Xsl^l~) versus Br(.B ^s'y) 
for RS parameter points with = 3 (0.5) and Mg(i) = lOTeV in black (green). More than 
70% (30%) of the model points he within the experimental 2a region. In fact, those numbers 
do not signihcantly change when switching off the electroweak dipole corrections Cf^{fib) 
and Cf^{yb)- The decay Br(R —)■ is mostly sensitive to the tree-level corrections 

in the RS model, justifying the discussion of the branching ratio in [13], where the 
corrections from the dipole coefficients have been neglected. Requiring that at least 10% of 
the RS parameter points he inside the 2a error margin yields the lower bound M^(i) > 3.3 TeV 
(6.2TeV) for =3 (0.5). 


5.3 Time-dependent CP asymmetry in 

In the SM, the left-handed structure of the weak interactions makes the emitted photon mainly 
left-handed in b decays (6 —)■ S'Jl) and right-handed in b decays (b —)■ s'yji), since the chirality- 
flipped Wilson coefficients are suppressed by ms/mb relative to the original ones. The helicity 
suppression of right-handed photons makes the time-dependent CP asymmetry dominated 
by R-meson mixing in the SM, irrespective of hadronic uncertainties. But in new-physics 
scenarios like the RS model there can be chirality flips on internal lines of a penguin diagram, 
such that the amplitude for a right-handed photon in b decays is no longer suppressed by 
mg/mb- Experimentally, the photon helicity can be accessed indirectly by using the time- 
dependent CP asymmetry in R° —)■ K*^'j decays, dehned as 


P(R°(t) ^ 77*%) - P(R°(t) ^ 77*%) 
T{B%t) 77*07) + P(R0(t) ^ 77*07) 


S'iftg, sin(Amst) — Cx^'y cos{AmBt ), 


( 82 ) 


where Ams is the mass difference between the heavier and the lighter neutral R-meson mass 
eigenstate. The mesons K*^ and K*^ are observed via their decay into the CP eigenstate 
Ksti^. The helicity suppression can be measured by SK*'y, which to leading order is given by 
[66, 67] 




K*-y 


|C'7g,(/if,)P -f \C7.y{fib)\‘^ 


Im 


e C7g,(/ife) C'7g,(/if,) 


(83) 


This observable is sensitive to the chirality-flipped Wilson coefficient Cj^i^yb)- The angle c/d is 
the phase of R° — mixing and has been measured in R —)■ J/tjj Ks decays to be sin^^ = 
0.682 ± 0.019 [54]. Due to the occurrence of Cj-yi^yb) in the numerator, the SM prediction for 
S'x*7 is suppressed by the ratio mg/mb and reads = (—2.3 ± 1.6)% [68]. The current 
experimental value = (—16 ± 22)% [54] still suffers from large uncertainties. 
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Figure 8: Corrections to Cg^io{fib) and Cg^io^fib) in the RS model for parameter sets 
with ?/* = 3 and Mg(i) = 10 TeV. The vertical dashed lines denote the lower and upper 
quartiles of the distributions, i.e. 50% of the RS points are included in the region in 
between the dashed lines. In the SM we take 4.07 and Cf^{^b) ~ —4.31. 

Similar plots are obtained for the imaginary parts of the Wilson coefficients. 


The right plot of Figure 7 shows the RS contributions to SK*'y and the branching ratio 
Br(.B Xsj) by the black (green) points for y* = 3 (0.5) and Mg(i) = 10 TeV. Gray and blue 
bands denote the experimental values with the la error margins for Sk*'-/ and Br(R Xg'y), 
respectively. Compared with the SM prediction the RS corrections can be significant due to 
the sensitivity of SK*'y on the imaginary part of which can receive large corrections 

in the RS model. On the other hand, the corrections are not signihcant when compared with 
the experimental result due to the large uncertainty. We observe that more than 95% (99%) 
of the points he within the experimental 2a region for = 3(0.5) and Mg(i) = 10 TeV. 
Requiring that at least 10% of the RS points lie within the experimental 2a regions of 
and Br(R —)■ Xg'j), we can derive the lower bound Mg{i) > 3.8TeV. 
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Figure 9: Shown is the direct CP asymmetry of i? —)■ as a function of the difference 
of the CP asymmetries of the charged and neutral B mesons. The light gray and 
blue bands denote the experimental la error margins and the area between the two 
horizontal dashed lines shows the la error margin of the SM prediction. The black 
(green) points represent RS points with y* = 3 (0.5) and Mg(i) = lOTeV. 


5.4 Direct CP asymmetry in B ^ Xs'j 

The direct CP asymmetry measures the difference between the rates of the decays B —)■ Xgj 
and B —)■ X^y and is defined via the ratio 

^ - Br(R ^ X,y) 

CP ^ > Bt{B ^ X,y) - Br(R ^ X,y) 

with a lower cut on the photon energy, which depends on the experiment. The experimental 
value = (1.5 ±2.0)% [54] is compatible with zero. Theoretically, the CP asymmetry 

is affected by perturbative “direct photon contributions”, in which the photon couples to a 
local operator mediating the weak decay in the effective low-energy theory, as well as non- 
perturbative “resolved photon contributions”, which account for the hadronic substructure of 
the photon. Taking both effects into account the SM prediction lies in the region —0.6% < 
— 2.28% [69] and is compatible with the experimental result. 

Let us now investigate the CP asymmetry in the RS model, where we can have additional 
contributions from chirality-flipped Wilson coefficients. The perturbative contribution [56] 

I ^ + Kz. i)] Im[(l + 

1^771 + X7l\ ^ 

- i Im|G8,C;, + + ^b{z, S) Im[(l + (85) 

+ ±i(2,i)|C2|Umle,]| 


(84) 
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is suppressed at leading order by a factor as{mb) arising from the strong-interaction phases. In 
the SM eg = Xu /At is the only source of a CP-violating weak phase, since all Wilson coefficients 
are real. In addition to this CKM suppression the SM result is further suppressed by the mass 
ratio 2: = ml/ml resulting from the GIM mechanism. The functions v{z), b{z,S) and b{z,S) 
can be found in [56, 69]. The non-perturbative contribution [69] 
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( 86 ) 


starts at leading order in Aqcd/''^?) and involves hadronic parameters with values in the range 
—330MeV < A27 < 525 MeV, —9MeV < A27 < llMeV and —17MeV < A78 < 190 MeV. For 
our analysis we choose the values A27 = 96 MeV, A27 = 1 MeV and A78 = 104 MeV. 

Adding both contributions (85) and (86) we obtain the black (green) RS points shown in 
the left plot of Figure 9, for y* = 3 (0.5) and M^(i) = lOTeV. We observe that all points lie 
inside the experimental 2a area. The effects decrease for smaller values of 7/*. In general, 
the corrections to the observable are too small in order to constrain the RS parameter 

space. 


5.5 The CP asymmetry difference in B Xs7 


Another observable is the difference of the CP asymmetries in charged and neutral R-meson 
decays. Its formula is given by [69] 


= Acp(5+ ^ X+7) - Acp(5° ^ X°7) 


A 


= 47r^Q;s(^^b) —- Im 


mb 


IC77P + IC77P J 


(87) 


where the hadronic parameter A78 is predicted to he in the range 12 MeV < A78 < 190 MeV. 
Recently the BarBar collaboration has published the first experimental result = (5 ± 

3.9stat 1.5syst)%, which is compatible with a null asymmetry difference at the level of la. For 
the analysis we take the average value A78 = 89 MeV. 

Figure 9 shows the predictions in the RS model for 7/* = 3 (0.5). For most of the points the 
corrections do not exceed the 1% level. All points are included in the experimental 2a error 
band. The effects are also decreasing for smaller values of 7/*. Currently, the RS parameter 
space cannot be constrained by the CP asymmetry difference 


5.6 Comparison with (almost) model-independent fits 

Several tensions at the 2-3a level in R —)■ K*MM angular observables have shown up in the 
collected data of the LHCb run during 2011 and 2012, including an integrated luminosity of 
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Figure 10: RS predictions for the real parts of the Wilson coefficients (left) 

and (right). Black (green) dots correspond to RS scenarios with = 3 (0.5) 

and Mg{i) = 10 TeV. The brown contours correspond to the la and 2a best £t regions 
obtained in [71] and [72, 78], respectively. The dotted lines represents the SM values. 


3 fb“^. One result is that the observable in the invariant lepton mass region 4.0 GeV < < 

8.0 GeV is only compatible with the SM prediction at the level of 3.7a [70]. Several model- 
independent theoretical analyses including also additional observables in B decays have been 
performed, showing that the deviations can be explained by new physics [71-77]. 

The full 4-body angular distribution of the decay B —)■ K*y'^y~ provides a sensitivity to 
the operators Q^-y, Qg, Qio and to the scalar and pseudo-scalar operators Qp, Qs and their 
chirality-flipped counterparts. In the RS model we neglect the corrections to 6 —)■ syPyT 
arising from the scalar operators Qp and Qs, since they follow from the tree-level exchange 
of the Higgs boson and are suppressed by the small lepton masses [13]. Here, we will focus 
on G 7 .^^g^io and their chirality-flipped counterparts. Unfortunately, in the literature there is no 
general £t where new physics is allowed to enter all Wilson coefficients at once and including 
both real and imaginary parts. Therefore, we have to consider certain scenarios where the 
underlying assumptions are partly violated in our model. This fact needs to be kept in mind. 

A general fit where new physics can modify Cjy^g^io{yh) and simultaneously 

has been performed in [71]. It includes experimental data on the B —)■ K*y~^y~ observables 
Pi, 2 , Pi, 5 ,6,8 and ApB in various bins, as well as data on Br(R —)■ X^y), Br(i? —)■ X^/i+p”), 
Bt{Bs —)■ y^y~), the isospin asymmetry Aj{B —)■ X^y) and SK*y Deviations of the Wilson 
coefficients from their SM values are dehned by G]^^(pfe) = Ci{yh) — Gf^(pfe). The best £t 
regions are given by (at 95% GL) 

C^^{yQ e [-0.06,0.01], C^^{yb) e [-1.8, -0.6], C^^iyQ e [-1.2,2.0], 

( 88 ) 

C^^{yQ e [-0.09,0.06], C^^{yb) e [-0.8,1.4], Gfo^(/ib) e [-1.0,0.8], 

showing that only G^^(/rb) is inconsistent with zero. The most economical scenario corre- 
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spends to a negative new-physics contribution to Cg{fj,b) with all other Wilson coefficients 
close to their SM values. Unfortunately, the RS model does not lead to such a large correction 
to Cg{fib), as is evident from the upper left plot in Figure 8. The corrections to the remaining 
Wilson coefficients are however compatible at 95% CL with the fit regions given in (88). 

In order to clarify the role played by some of the Wilson coefficients some constrained 
scenarios have been considered, in which only two Wilson coefficients at a time are assumed 
to receive contributions from new physics. While this assumption does not hold true for the 
RS model, we still like to compare our results with the two scenarios where only Cj^i^Hb) and 
Cgi^fib) or Cg^jjib) and Cio(/ife) are assumed to be modihed by (real) new physics contributions. 
We do not consider modiheations of Cg{fJ,b) and Cio{fib), since the corresponding corrections 
in the RS model are very small. In the left plot in Figure 10 we consider the 
plane. The brown contours correspond to the la and 2a best £t regions obtained from the 
analysis in [71]. The best fit points are given by Cf^{fib) ~ —0.02 and Cg^{fib) ~ —1.5. 
We observe that almost none of the RS points touch the best fit region, since Cg{fib) does 
not receive large enough corrections. As a second scenario, we consider in the right plot 
of Figure 10 the plane, where the best fit regions are obtained from a global fit 

performed in [72, 78]. The fit includes 88 measurements of 76 different observables, including 
B —)■ angular observables and branching ratios as well as the branching ratios of 

B —)■ B —)■ Bs —)■ B —)■ K*'y, B —)■ and Bg —)■ We 

observe an anti-correlation between C^^{fj,b) and (/i;,), which is also favored by the best 
fit regions. While the corrections (/i;,) can be quite large for a few points in parameter 
space, the corrections C^^{fj,b) are too small to reach the best fit regions. 

6 Conclusions 

We have investigated the electro- and chromomagnetic (quark) dipole coefficients for b ^ s'j 
and b ^ sg transitions in the minimal Randall-Sundrum model with a brane-localized Higgs 
sector. We have derived integral expressions for all contributions arising at one-loop order using 
5D fermion and gauge-boson propagators and retaining the full dependence on the Yukawa 
interactions. The expressions are formally valid to all orders in in contrast to [16], 

where the Yukawa interactions were treated as small perturbations. Our final results involve 
one momentum and two extra-dimensional integrations, and each integrand contains one 5D 
gauge-boson and one 5D fermion propagator. 

By analyzing the UV behavior of the 5D propagators we have confirmed the finiteness of 
the penguin loops, as shown first in [16]. In addition, we have derived expressions in the KK- 
decomposed (4D) theory and shown analytically and numerically that the dipole coefficients 
coincide in both pictures, presenting a highly non-trivial cross-check of our calculations. We 
have derived approximate formulas for the KK contributions to and and have 

shown that the dominant corrections originate from the lU^-boson penguin diagrams. More 
precisely, when working in Feynman-’t Hooft gauge the dominant corrections stem from the 
parts of the diagrams which involve the scalar component of the 5D gauge-boson held and the 
charged Goldstone bosons from the Higgs sector. We find that for not too small values of the 
anarchic (5D) Yukawa matrix entries the latter contributions dominate and the size of the KK 
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corrections to the dipole coefficients increases proportional to yl- In contrast to [17], we have 
not found a signihcant contribution of the triple gluon vertex penguin diagram on C^g and C^g. 
In agreement with [17], we have observed the general pattern that the chirality-flipped Wilson 
coefficients CY^y^sg receive larger corrections than since the left-handed bottom-quark 

prohle is more localized towards the IR brane. For the Higgs and Z-boson penguin diagrams 
we have obtained results in the brane-localized Higgs and narrow bulk-Higgs scenarios. For 
I/* > 1 both contributions cancel to good approximation for the case of a brane-localized Higgs, 
while there remains a non-zero (zero-mode) contribution for the case of a narrow bulk-Higgs. 
A similar observation was reported in [20, 48-50] for the case of the leptonic dipole coefficients. 

In our phenomenological analysis we have RG evolved the dipole coefficients to the R-meson 
scale Hb- We have then performed a numerical scan of the RS parameter space with anarchic 5D 
Yukawa matrices and investigated the branching ratio Br(R —)■ Xg'y), the time-dependent CP 
asymmetry S'x* 7 , the direct CP asymmetry and the CP asymmetry difference 

all of which are sensitive to corrections to the dipole coefficients. Currently, the observables 
Br(R —)■ Xg'y) and can be used to constrain the RS parameter space. Requiring that at 
least 10% of the RS model points lie in the 2(7 experimental error margins, we can derive a 
lower bound on the mass of the hrst KK gluon resonance of M^(i) > 3.8 TeV for Yukawa matrix 
entries bounded from above by y* = 3. For smaller values of the bound gets weaker. We 
further discussed the branching ratio Br(R —)■ Xsl~^l~), which is dominated by the tree-level 
corrections to C'io(/ife) in the RS model. For this observable we can derive a lower bound of 
Mg(l) > 3.3 TeV (6.2 TeV) for y* = 3 (0.5), showing that the RS corrections increase for smaller 
values of y^,. Finally, we have compared the Wilson coefficients and (^ 77 , 9,10 with the 

results of (almost) model-independent hts performed in [71, 72, 78]. In general, the tree-level 
corrections to Cgi^yb) are too small in the RS model in order to cover the best £t regions. 

Node added: While this paper was under review the work [79] appeared, which contains a 
detailed analysis of the decay B —)■ in the minimal and custodial RS model with an 
IR-localized bulk Higgs. Their implementation of the scalar sector includes contributions from 
the Higgs zero mode and its KK excitations. In contrast to our approach the Yukawa inter¬ 
actions are treated as perturbations and results for the one-loop penguin diagrams to leading 
order in are derived. In addition to our work, the mixing of the dipole operators 

with four-quark operators obtained by integrating out KK gluons and the Higgses are in¬ 
cluded. Comparing our results for the branching ratio Br(R —)■ X^y) in the minimal RS model 
with the analysis of [79], we find that they are of similar size. 
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A Summary of Feynman Rules 

All particle momenta are flowing into the vertex and amplitudes are denoted by A. 


5D Theory 

We begin with the vertices that couple two quarks and one boson. Each vertex is accomplished 
additionally by an integral dt and we obtain for the corresponding amplitudes (q = u, d) 


Cw \/1 > 


A{qA5q,qaG5qb,qZ5q,qW^q'} _ j rpc 95 Ql 95 \ 


v^At)PL + v^At)PR , (A.i) 


= _± _ 1 ) + Af 


where P 12 + P 21 , 9 l = — Qdsl,, e = is the 4D electromagnetic and 

9s = 9s 5 / is the QCD 4D gauge-coupling. In the hrst two lines the subscript B of Pb, 
(t) and Vg±(f) must be replaced by the corresponding boson label B = A,G, Z,W on the 
left side®. The t-dependent functions are given by 


VA^Psit) - ^ 1 > 


^w+it) - ^ 




V.M = ^ 




Lm^ 


*'ud 


Pw + 




Lm^ 




du 


Pz+{1 


9l\^ 

9l) Lm% 


2 

KK 




- -VA^pt) , 

u,,-*(«) = -V„p) 

W(«) = -v„p) 

VzM = -VzM 








(A.2) 

where gl^ = —Qqsl^ and = Y^ P 12 — Y^^P 21 . We continue with the Feynman rules 

for the triple gauge-boson vertices, where we assume that the external photon line is a zero 
mode. The corresponding amplitudes read 




W^{p)APq)W^ik) ^ ie5^ [, 

^ Lt 


- py + g^yp - k)^ + g^^{k - qf 


jW^{p)APq)Wf{k) ^ ^g/g ^TTMKK 

L 


dY^ +5yt-i) 


p ’ 




2ti 




e 

E ’ 


^The superscripts of are only relevant ior B = W and can be ignored otherwise. 
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^G^(p)GS(g)G^(fc) = ^.,5^ "27r l ^ ^ 

^G^(p)GS(<?)G'g(fc) = I jabc ^TrMKK ^Gg ^ 

^G|(p)GSGg(fc) ^ ^afec (p _ ^ 

Zy t 


(A.3) 


where the superscript of the t-derivatives indicates the held it should act on. 


4D Theory 

In the KK-decomposed (4D) theory the amplitudes for the vertices coupling two quarks and 
one boson can be summarized by 


j\^{'lAij,q,qaG'^^qb,qZfj,q,qW;tq'} _ 
_^{q‘PAq,qa‘Paqb,wzq,Wwq'} = 

= -i 


(■ rr^c ^95 91 Wb} 7 a * 


^nmk Pl + Vumk Pr 


Cyo V2^ 

{gDnk Pl + {gDnk Pr 


y'PB p , yVB p 
n.m.k. ' n.rti.k. 


(A.4) 


where n, m, k are the mode-numbers of the anti-quarks, bosons and quarks respectively. The 
labels B^ifs on the right side must be replaced by the corresponding boson (label) on the left 
side. The vector and scalar overlap integrals for the neutral gauge bosons are given by 


r = 7^ 
O = 72^ 


it xi(t) Qf\t) Pb S?V) . . 

it Sg”’'(*) VbM Qt\t) . 7„« = VZ^k„R.v,,-^v,. 


(A.5) 


where B = A,G,Z. The structures VBg(t) are dehned in (A.2). For the boson we hnd 
the following overlap integrals 


v,H = 7S 


VnZ = 7 S 


= 72^ 


<itx;:(t)«F'(t)Ui-i>r(‘)- vzi = vzi\b»v, vzi = vii\L»R, 
ufUt) w,+(() vf{t) . zL = uAl 


■)(^) / 


rW+i 


rW± 


xW±i 


Loi?,VV+^V^+ ’ 


(A.6) 


dt 


-ktdtx'^jt) (n)t 




mw„ 


yVw _ yV’w I 

^ nmk ^ nmk\ LtAfR^V -—rV - 

Wc Wc 


The Goldstone-boson contributions from the Higgs sector are contained in the brane-localized 
terms of Vg±(f) in (A.2). We can simplify the scalar overlap integrals by performing a partial 
t-integration noting that boundary terms at t = e, 1 vanish, since we work with the regularised 
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5-function (19). We can apply the equation of motions for the quark prohles such that the 
terms with the 5-functions cancel. For instance, in case of we can use that 


dt 


4”'t(i)p+cr(i) 


(fc)/ 


m 


F+Cf («) - PPit) 

Mkk 


(A.7) 


where the last term in (A.7) cancels the 5-function appearing in Vy^+{t). Repeating the steps 
for the remaining cases we hnd 


yVB ^ yB _ yB y^B ^ yB _ ^qk t/B 

^ nmk ^ ^ nmk ^ ^ nmk ’ ^ nmk ^ ^ nmk 










{B = A,G,Z) 


yVw ^ yW+ _ ^dk yW+ yVw ^ yW+ _ ^dk yW+ 

^ nmk ^ 'nmk ^ ^nmk ) ^ nmk ^ ^nmk ^ ^ nmk i 

mw^ mw^ mwm 


(A. 8 ) 


yVw _ ^d„ T^VC- _ ^Uk yW- y'Pw _ ^d„ yW- _ ^^k yW- 

^ nmk ^ 'nmk ^ ^ nmk ) ^ nmk ^ ' nmk ^ ' nmk ) 

mw^ 


where all scalar overlap integrals can be expressed in terms of the vector overlap integrals. 
Concerning the quark-Higgs-quark couplings in (A.4) the overlap integrals are given by (g = 
M, d) 


^A)nk = -^j^ 1) Si’'(<) SS’b). Wnk ^ [(gl)tn]'■ 


(A.9) 


In the brane-localized Higgs scenario with = Y^ we can perform the f-integration analyt¬ 
ically and find 


^ oWtn-'i " 


Mu = ^sr(i-)p.-,2x, 


us®(i"). 


(A.IO) 


where ^ Q{YiyyP, Yq = (taiih Xg/Xg)l^ and {gl)nk = [(9a)rn]^- We eontinue with the 
triple gauge-boson vertices where photon or gluon zero modes are attached. We obtain 
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Trr L 
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— py — ky -t- ? 7 '*'"(/;; — q)^ 


^nk ■ 
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ley, 

V2717 


^nk 
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(A.11) 


^nk 
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^TTlQ^Yik 1 
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^rG,l(p)G'«°Ug)rG,l(*:) = (p - ky ■ 
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B Solutions for the 5D quark propagator 


For details on the procedure of calculating the 5D quark propagator (12) in the mixed position- 
momentum space we refer the reader to [29], where the solutions for the propagator functions 
F; A;|;) and F; A;|;) have been derived (with Euclidean momentum k% = —k^). 

Here we extend their results and also include the solutions for F; k%). The solution for 

k%) can be obtained by complex conjugation = [A?ji(f,t';A;|)]^ We 

begin with the solutions in the minimal RS model for a brane-localized Higgs sector where 
the regulator 77 of the regularized (5-function in (19) fulhlls the constraint 77 T/^n/Axev, 
implying that rjkE y*v. As a consequence the 77 dependence drops out of the solutions for 
the propagator functions. In this limit rjkE ^ y*v we obtain for q = u,d: 


^LL 
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^LL 


^LL 


Ag ,22 

^LL 
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D?{kE,t') D'^{kE,t 
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-yy \ D^{kEy Zg 1 DfCkE,t') 
kEMKK ID^ikE, l)^ + Zq Rq D^CkE, 1) 


+ 


D^{kE,ty 

D^{kE,l) 


L^{kE, l,t>) 


5 


48 



q,12 _ D 2 i^E,t) Zg 1_1 Dl{kE,t’) 

~ kEMKK D^{kE, l)^ + Zq Rq gYg^ Df{kE, 1) ’ 

q,2l _ -Vt^ Dl{kE,t) 1 Zq 1 D^{kE,t') 

/cbMkk Dl{kE, 1) gYq 1 + Rq D^ikE, 1) ’ 


^ 9,22 

^RR 


-VtF 


Dl{kE,t) 1 


DlCkE,t') DlCkE,t^, 


Ll{kE, 1 , 


kEMKK [DlCkE, 1 ) g% l + ZqRQ gYq^ DlCkE, 1 ) DfikE, 1 ) 

For the sake of readability we have suppressed the arguments of the propagator functions 
F;/c|;) for A, B E {L,R} and of Zq{k%), Rq^qikE) dehned in (29), (30). The mod- 
ihed Yukawa matrix is dehned by = (tanhXg/Xg) Yq with Xq = QiYqYkf/^ and q = 

u/(\/2Mkk)- We have used the abbreviations kE = /cb/Mkk, t> = Max(t, F) and t< = 
Min(t,t'). In (B.l) we used the functions Df’'^{kE,t) [29] which are related to the more gen¬ 
eral functions Df’^{kE,t,t') dehned below via Df’‘^{kE-,t) = Df’^{kE-,t,e). The generalized 
functions read 

D-^/kE,t,t') = I_,.{kEt')h 


.{kEt) 




l(M')/_c.±l(M), 


D^^^{kEAA') = I-aA + ^(^Et') Ic^^AkEt) - Ici{kEt') I_^^^i{kEt) , 
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2 COS { eca ) 


DfikEAA')] i = l,2,3,4, 


(B.2) 


(B.3) 


for A = Q,q. We introduced the shorthand notation Ijj = l — rj and the bulk-mass parameters 
are denoted by CQ g. 

Next we focus on the case of a narrow bulk-Higgs where the regulator r] takes values in the 
range y^v/ATev ?7 I/*u/Mkk, implying that the propagator solutions explicitly depend 
on the product rjkE- For the calculation of the scalar contributions in Section 3 we need 
to evaluate the 5D propagator functions in the region near the IR brane, where the extra- 
dimensional coordinates take values in the range t,t' G [ 1 , 7 , 1 ]. Focusing on the solutions in 
this region we obtain the results (t,t' G [ 1 , 7 , 1 ]): 
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Again we have suppressed the arguments of the functions for the sake of readability. The 
step functions is denoted by 9(t — t'). The t-dependent functions 5(t), C{t) are defined in 
(23) where we have to replace Xq by Sq as dehned in (25), analogously for S{t) and C{t). In 
(B.4) the modified Yukawa matrix is dehned hy Yq = (tanh Sq/Sq)Yq. We also introduced 
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new structures (g = m, d) [29] 


Nf(k%) = l + Zf'{kl)+tikE 
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Nfikl) = 1 + Zfikl) + ,,kE 


Zf(kl)Rg'Ck 


coth Sa tanh Sa 

I—^ ^ 




(B.5) 


with 


Zfikl) = S^^Y, R,ikE) Y^RqikE) . Zf(k:i) = fY, R,(kE) Y^ ^ Rq^e) 


(B.6) 


C Ultra-violet behavior of 5D propagators 

This section discusses the behavior of the 5D propagator functions in the brane-localized Higgs 
scenario for large momenta /ce 3> Mkk/^ exceeding the effective Planck scale at each point in 
the extra dimension. The results are used to show the hniteness of the penguin diagrams and 
to calculate the boundary terms in Section 3. 


Gauge-boson propagator functions 


Expanding the scalar and vector parts of the gauge-boson propagator functions in Euclidean 
momentum space we hud to leading order (/ce 3> 
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(C.l) 


for (subscript) B = H, G, W, Z. In case of the massive gauge bosons B = W,Z the scalar 
propagator function includes the contributions from the hfth component of the 5D gauge- 
boson helds and from the corresponding Goldstone-bosons in the Higgs sector, which gives 
rise to the term proportional to Lm‘^/kEM-K_K- This term is absent in case of the photon or 
gluon scalar propagator {B = H, G). Integrating (C.l) with a (well-behaved) function 
along both extra dimensional coordinates we can show that 


dtdt' Bsit, t'] A;|) /(f, t') 


dtdt' B~{t,t']kl) 


ATikl J, 


dtf{t,t), 


CkE > 1/e) (C.2) 


87re2/c| 


dtf{t,t). 


We cannot prove (C.2) in general but we checked analytically that the relations are valid for 
the functions relevant in the calculations of this paper. Relations (C.2) imply that for large 
Euclidean momenta the propagator functions can be effectively replaced by the ^-function 
6{t — t') apart from a constant factor. 
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Quark propagator functions 

The results for the 5D quark propagator functions in the brane-localized Higgs scenario are 
listed in (B.l). The solutions contain several functions that have a simple form when expanded 
for large Euclidean momenta. To leading order we hnd the expressions {kE 3> 1/t) 


D^,ACkE,t) 


Vi 

Lt^2{kE, 1, t) = -Lf^^{kE, 1, t) 


I _|- g2fcB(e-t) 


DUkE,t) 


’ D^^,{kE,l) Vi 

gfcB(i-t) 


I _ g2fc£;(e-t) 


, RaVe) ~ 1, 


2Vi 






(C.3) 


for A = Q,q and q = u,d. The expansions are independent of the bulk-mass parameters. Using 
the above expressions in case of the propagator function f; k%) we hnd {kE > 1/t, 1/t') 


'^''l ^E 
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^eMkk 1 + g^YY^ ^ 


^ifVt'VE. 




-kE{2-t-t') 

-- 


-q-^q 
1 


X q- g2fcB(e-t)1 Tq _ Q2.kE{e-t') 


(C.4) 


/ceMkk ^ 1 + g^YqY;! L 


X _ g2fc£(e-i)l ["x q- g2fcis(e-i') 


^-kE{t>—t<) 

2/c_bMkk 


X — g2fc£(«-*<) 


1 _ 1 - 
1 + 


at leading order in k^ . For a (well-behaved) function f{t,t') we can show that 


dtdt'A^^(t,t';/c|)/(t,t') 


kl Je 


dtf{t,t), (fc£;>l/e) 


(C.5) 


implying that the propagator function behaves like the ^-function 6{t — t') for large momenta. 
For functions that are localized near the IR brane the equation is already a good approximation 
for kE 3> 1. Equation (C.5) is also valid in case of the propagator function ^RRVt'Vl)- We 
continue with the chirality-changing propagator function \\^{tA']k\). Using the coupled 
differential equation [29] 


^ j ^e) ~ "14kk( df + ■M.q{t )) A|^^(t, t ; kE ), 

we can show for a (well-behaved) function f{t,t') that 


(C.6) 


dtdt’ Alj^{t,t'-,kl) f{t,t') = Mkk / dt - A|^(t, 1; A;|)/(t, 1) + A^^(t, e; A;|)/(t, e) 


+ Mkk / dtdt'Al^{t,t';kl) d, +M,{t') f{t,t’), (C.7) 
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where we have performed a partial integration in t'. In the large Euclidean momentum region 
we can show that {kE 3> 1/e) 




, P_ 
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-q 


12 


-qY, 
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21 
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(C.8) 


dt e; k%) fit, e) —-y P+ fie, e), 


k% 


where higher order terms are suppressed at least by k^^. The hrst relation is approximately 
valid already for /c^; 3> 1 if fit, 1) has most of its support near the IR brane. Using (C.8) in 
(C.7) we hnally hud for a (well-behaved) function /(t, T) {kE 3> 1/e) 


dtdt’ Al^{t,t’-,kl) fit, t’) 
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kl 
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(C.9) 


where f'it,t^) = lim^^tio /(t, s) is understood as a limiting procedure. An analogous 
equation can be derived for the propagator function A^^(t, T; kff). 


D Loop functions 

In the KK-decomposed 4D theory the Wilson coefficients in (57) involve the loop functions 
h^iix) and iQ^uix). They are dehned by the integral representations 
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(D.l) 


where = ml/ml- Note the different arguments of the loop functions /a, 4 , 6 , 7 ( 2 ^) and /8_ii(a;). 
Performing the momentum integrals the loop functions explicitly read 
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